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A novel method of studying e+e~ annihilation into hadrons using initial state radiation at e+e" 
colliders is described. After brief history of the method, its theoretical foundations are considered. 
Numerous experiments in which exclusive cross sections of e+e~ annihilation into hadrons below 
the center-of-mass energy of 5 GeV have been measured are presented. Some applications of the 
results obtained to fundamental tests of the Standard Model are listed. 
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FIG. 1 The lowest-order Feynman diagram describing the 
process of e^e~ annihilation into hadrons. 



FIG. 2 The lowest-order Feynman diagram describing the 
initial state radiation process e^e~ — >■ 7-I- hadrons. 



I. INTRODUCTION 



A. Why is low energy e+e annihilation interesting? 

Studies of low energy e^e~ annihilation into hadrons 
are of great interest for theory and have numerous appli- 
cations. According to current concepts, e^e^ annihila- 
tion into hadrons proceeds via an intermediate virtual 
photon which produces a pair of quarks, qq, followed 
by the hadronization of quarks into observed hadrons. 
This process is described by the lowest-order Feynman 
diagram shown in Fig. [T] When the initial energy of 
e^e~ , or equivalently of the intermediate virtual photon, 
is large enough, the process of hadronization is well de- 
scribed by Quantum Chromodynamics (QCD). At small 
energies, lower than 2-3 GeV, produced hadrons are rela- 
tively soft and intensively interact with each other form- 
ing hadronic resonances. At the moment QCD fails to 
describe this energy region. Because of that, it is vitally 
important to gain sufficient information from experiment 
to be used as an input to various QCD-based theoretical 
models. QCD sum rules are an example of how measure- 
ments of total and exclusive cross sections can be used 
to extract such fundamental parameters of theory as the 
stron g coupling constant a s, quark and gluon conden- 
sates ()Shifman ali\2QQ^ . 



Precise knowledge of vacuum polarization effects 
based on the total cross section of e+e" annihi- 
lation into hadrons is necessary to estimate the 
hadronic contributions to the running fine-structure 
constant and thus determine its value at the Z bo- 
son mass, a{M'^), a key compon ent of the high- 
precision tests of t he St a ndard Model (lActis et al\ . | 2010|: 
Burkhardt et al!, '1989; 'Burkhardt and Pietrzyk", '200I 
Eidelnian and Jcgcrlchncr, 1995; Hagiwara et at, 200^. 

Improvement of the precision with which the to- 
tal cross section of e+e~ annihilation into hadrons 
is known is also needed for a more accurate es- 
timation of the hadronic contribution to the muon 
moment since it is one of 
in a search for New 



IS 

of the 
anomalous magnetic 
the cruci al limiting fa ct ors 



There is an important relation between spec- 
tral functions in e+e^ annihilation into hadrons 
with isospin / — 1 and corresponding t lep- 
ton decays based on conservati on of vector current 
(CVC ) and isosp in symmetry ( Thacker and Sakurai 
I1971I: iTsai ll97lD . While first detailed tests of 
such relations showed satisfactory agreement between 
such spectral functions jEidelman and Ivanchenkol[l99ll : 
[Kawamoto and Sandal . 1 19781) . higher accuracv reached in 
both e+e~ and t lepton sectors revealed possible system- 
atic effe cts not accounted for in the e+e" and/or r exper- 



iments ( Davier et aU [2003al lbh. Understanding of these 
effects is crucial for improving the accuracy with which 
the hadronic contributions to the muon anomalous mag- 
netic moment can be estimated from r decays to two 
and f our pions as was first suggested by ( Alemanv et al\ . 
I1998D . 

Detailed measurements of the energy dependence 
of various exclusive cross sections allow to improve 
our knowledge of vect or mesons and l ook for new 
states, both of ligh t (iDruzhininl . l2007l ) and heavy 
quarks (|Eichten et all l2008l) . 



B. Idea of ISR 

In e^e~ collider experiments exclusive and total 
hadronic cross sections are usually measured by scanning 
the accessible energy range. The process of e^e~ an- 
nihilation is accompanied by emission of one or several 
photons from the initial state. The lowest-order Feyn- 
man diagram describing initial-state radiation (ISR) is 
shown in Fig. [2] The quantity measured directly in the 
experiment is the visible cross section 



N 



(1) 



Phvsi cs (iBennett et all. 120061: 



119611 : iGourdin and de Rafael) . 119691) 



Bouchiat and Michel) . 



where N is the number of selected events of the pro- 
cess e^e^ — > hadrons + 717, n = 0,1,2,..., and L is 
the integrated luminosity of the collider collected at the 
center-of-mass (cm.) e+e~ energy 2i?o. The visible cross 
section can be related to the Born cross section ctq cor- 
responding to the lowest-order diagram of Fig. [1] via the 



3 



integral ( Kuraev and Fadiiil . Il985l) . providing the 10 
accuracy: 



e{s,x)W{s,x)aa{s{l-x))dx, (2) 



where s = AEq , x is an effective fraction of the beam en- 
ergy Eq carried by photons emitted from the initial state, 
™min is the minimal possible invariant mass of the final 
hadrons, e{s,x) is the detection efficiency for the process 
e'^e^ — > hadrons + as a function of x and s. The 
so-called radiator function W{s,x) taking into account 
higher-order QED contributions, in particul ar, from the 
diagr am in Fig.[2j is fully calculable in QED ( Actis et~ai] . 
I2OIOD . Due to the photon emission from the initial state 
the visible cross section depends on the Born cross sec- 
tion at all energies below the nominal e'^e~ cm. energy 



2Eo. 

In conventional scanning experiments the influence of 
ISR is suppressed by the requirements of the energy and 
momentum balance between the final hadrons and the 
initial e+e^ state. In this case the detection efficiency has 
X dependence close to the step function: e{s,x) — eo{s) 
for x < xq, and zero for x > xq. At small xq , the equation 
([2]) can be rewritten: 



= eo(s)ao(s)(l + (5(s)). 



(3) 



where 1 + (5(s) is the radiative correction factor, which 
takes into account higher-order QED corrections. To cal- 
culate this factor it is necessary to know s dependence of 
(To in the range from s(l — xq) to s. For slowly varying 
cross sections, d is about 10%, and can be determined 
with accuracy better than 1% using existing data on the 
cross section energy dependence. Thus, in scanning ex- 
periments, from the data collected at the cm. energy 
the cross section ao{s) is determined directly. 

Another approach is also possible. Equation Q can 
be rewritten in the differential form: 



d(Tvis(s,™) 2m 
dm s 



e(s, ■m)W{s, x)CTo(m), 



(4) 



where we have made a transformation to the variable 
m — ^ys(l — x), the invariant mass of the hadronic sys- 
tem. At non-zero x the dominant contribution to the visi- 
ble cross section comes from the one-photon ISR (Fig. [2). 
With the inclusion of the ISR photon momentum into 
the selection conditions on the energy and momentum 
balance, the non-zero detection efficiency for ISR events 
can be obtained in a wide range of the hadronic invariant 
mass. So, from the measurement of the mass spectrum 
for the process e+e~ — ^ hadrons -I- 7 at fixed cm. en- 
ergy ^/s the cross section ao{m) can be extracted in the 
invariant mass range from threshold to the mass close to 

The idea of utilizing initial-state radiation from a 
high-mass state to explore electron-positron processes 
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FIG. 3 The relative probability for the ISR photon to be 
emitted into the polar angle range 60 < 9 < 180° — 9o for 
three representative values of x. 



at ail energi es below tnat s t ate was outlmea long 
ago m Refs. (jBaier and Fadinl . Il968l: H aier and Khoze l. 
Il965h . A possibility of exploiting such processes 
at high l uminosity (j)- and i3 - factories was discussed 
in Refs. (lArbuzov aL . 1998 : Benavoun aLl. 19991 : 



iBinner e/aLl . ll999t iKonchatnii and Merenkovlligga^ and 

motivated studies described in this paper. 

Analysis of ISR events at e+e^-factories provides inde- 
pendent and contiguous measurements of hadronic cross 
sections in the low-energy region and also contributes to 
the spectroscopy of low-mass resonances. 



C. Calculation of ISR and accuracy 

In the lowest order (Fig. [5]) the probability of the 
initial-state radiation of the photon wi th the energy xEn 
and t he polar angle 9 is as foll ows (jBaier and Khoze , 
119651 : iBonneau and Martini . 1 197lh : 



wo{d,x) 



a 

TTX 



(1 



■ sin^ 9 



sin^ 9 



■ cos^ 9 



Ami (1 ~ 2a;) sin^ 



(smH+^cos^9 



(5) 



where a is the fine-structure constant, and me is the 
electron mass. 

The ISR photon is predominantly emitted at small an- 
gles with respect to the beam axis. In Fig. |3]we present 
the dependence of the function Wo(9o, x)/Wo{0, x) on the 
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FIG. 4 The relative probability for the ISR photon to be emit- 
ted into the polar angle range 30° < 9 < 150° as a function 
of the e^e~ cm. energy for three representative values of x. 



polar angle limit where 



Wo(6'o,x)= / wo{9,x)sm9de. 



(6) 



The integration is performed for three values of x at 
2Eo — 10.58 GeV, the cm. energy of B-factories. It 
can be seen that the angular distribution of the ISR pho- 
ton weakly depends on x and that a considerable frac- 
tion of the photons is emitted at large angles. In the 
next section we will discuss two approaches to study ISR 
events, a tagged and untagged one. In the tagged ap- 
proach the ISR photon should be detected, i.e., emitted 
at a large angle, into the fiducial volume of the detector. 
At B-factories {2Eo = 10.58 GeV) about 10% of high- 
energy ISR photons have 30° < 9 < 150°. This angular 
range approximately corresponds to the fiducial volume 
of the electromagnetic calorimeter of the BABAR detec- 
tor. The fraction of the large-angle ISR increases with de- 
crease of the energy as shown in Fig. 21 The compact ex- 
pressions for Wq can be written for two practically appli- 
cable cases. For the range of integration 9o < 9 < tt — 9o, 
9o > me/^/s 



Wo{9o.x) = — 

TTX 



{2-2x + x'^)\n 



1 + cos 6*0 
1 — cos 9o 



-X cos 00 



(7) 



For the full range of polar angles Q < 9 < tt 



Wo{Q,x) = —{\n^-l){2-2x + x'). (8) 



TTX mt 



The formulae given above describe ISR processes in the 
lowest QED order. To estimate a contribution of higher- 



FIG. 5 The mass (m = 2Eo\'T—x) dependence of the rel- 
ative difference_betw^n tlw rad^ function W{x) from 
Ref. (H uraev and Fadinl . Il985l ) and the lowest-order function 
Wo{0,x) for 2Eo = 1.02 GeV. 



order diagrams (loops and related to extra photon emis- 
sion) the function W{x) from Ref. ( Kuraev and Fadinl . 
119851) can be used, which takes into account soft mul- 
tiphoton emission and terms in the leading logarith- 
mic approximation. In this approximation the accuracy 
AW/W is expected to be better than 1%. The rela- 
tive difference between W{x) and Wo{0,x) as a func- 
tion of the invariant mass of the final hadronic system 
is shown in Fig. [S] for 2Eo — 1.02 GeV, the cm. energy 
of the ^-factory in Frascati. It is seen that the radiative 
correction to the lowest-order radiator function reaches 
15%. It should be noted that the size of the radiation 
correction depe nds on experimental conditions. For ex- 
ample, in Ref. ( Aubert et aLl . [2C)06a[) the function W{x) 
is calculated at 2Eo = 10.58 GeV with conditions that 
the highest-energy ISR photon has a polar angle in the 
range 20° < 9 < 160° and that the invariant mass of the 
hadronic system combined with the ISR photon is greater 
than 8 GeV/c^. The latter condition restricts the max- 
imum energy of extra photons emitted from the initial 
state. With these conditions the radiative correction fac- 
tor 1 + ^ = W{20°,x)/Wo{20°,x) is close to unity with 
the maximum deviation d of about 2%. 

To provide accuracy better than 1% required for 
the measurement of the exclusive hadronic cross sec- 
tions at low energies, the calculation of the radiator 
function should include the higher-order radiative cor- 
rection, in particular, due to emission of extra pho- 
tons. Several theoretical papers are devoted to study 
radiative corre cl ^ions to ISR proce s ses, f or example , 
(lArbuzov et all Il998l: iBinner et all Il999t ICzvz et all 



120031: iKhoze et all l200ll 120021: iRodrigo et all 120021) . 

The a pproaches of R efs. ( Binner et all 1999HCzvz et all 
120031 : IRodrigo et all \20m allow one to develop genera- 
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tors of Monte Carlo (MC) events a nd are used in anal - 
yses of experimental data. In Ref. ( Binner et al\ . Il999l ) 
the photon emission at large angles only is considered; 
radiative corrections are calculated in the leading loga- 
rithmic a££roximation_with^^ structure function tech- 
nique (|Caffo et all ll994L Il997t ). The accuracy of the 
method is determined by neglecting sub-leading con- 
tributions and estimate d in Ref. (R odrigo et al., 2001) to 
be ab out 1%. In Refs. (ICzvz et all 120031 IHodrigo lit ali 



I2002D the one-loop corrections and exact matrix element 
for emission of two hard photons are calculated. The ac- 
curacy of this next-to-leadin g order (NLO) calcul ation is 
estimated to be about 0.5% ( Rodrigo et aLl . [2002l ) due to 
the higher-order effects. 



D. Monte Carlo generators 

The calculation of the radiator function is usually per- 
formed by the Monte Carlo method. A special com- 
puter code referred to as an "event generator" provides 
events (sets of the four-momenta of the final particles) 
distributed over the phase space according to the matrix 
element squared of the process under study. The phase 
space can be restricted by some conditions on the an- 
gles and energies of the generated ISR photons. These 
conditions should be looser than the actual experimental 
conditions used for event selection. 

The interaction of the generated particles with the de- 
tector and the detector response are then simulated. In 
modern exper iments the detector sim ulation is based on 
the GEANT4 (|Agostinelh et aLl . l2003f ) package. The sim- 
ulated events are reconstructed with the program chain 
used for experimental data. The detection efficiency is 
determined as the ratio of the mass spectrum of simu- 
lated events that passed selection criteria to the spectrum 
of generated events. 

Most of ISR analyses discussed in this paper are based 
on two event generators. Historically, EVA was the first 
ISR Monte Carlo generator. The AfkQed package used 
in the BABAR experiment at the S LAC i3 -factory is a 
development of the E VA generator (| Binne r et ai, 1999; 
ICzvz and Kiihnl . l200l[ ) initially designed to simulate ISR 
production of the 2tt and 47r final states with an ISR 
photon emitted at large angles. The soft-photon radia- 
tion from the initial st ate is gener a ted w i th th e struc- 
ture function method (jCaffo et all Il994 Il997t ). Two 
extra photons are emitted in the directions of the ini- 
tial electron and positron. The program has a modu- 
lar structure allowing to implement easily new hadronic 
modes. The AfkQed package includes generation of 
27r, Stt, 47r, 57r, Gn, 7y7r+7r~ states, modes with kaons 
KK + mr n = 0, 1, 2, 3, 4, and protons pp, pp2iT. The 



generation of the process e 



7 is also included 



into the AfkQed package. For this process both initial- 
and final-state radiation (FSR) diagrams and their inter- 
ference are taken into account. For the charged particles 
the final-state radiation is generated using the PHOTOS 



package ( Barberio et a/.l . [l99lf) . 

The Phokhara event generator is used in the BABAR 
and Belle experiments at the S-factories, and in the 
KLOE ex periment at the ( ^-fact ory. Its latest version 
6.1 (PHOKHARA web sitd . [2009l )_ includes generation of 
the 277, 37r, 47r, KK, pp, and AA hadronic states, and 
the process e"'"e~ — )■ i X^jiT^. The initial-state radiation 
is ge nerated in NLO ( Czvz et~al\ . [200 3: Rodrigo et ad 
|2002| ). i.e., one or two photons can be emitted by the 
initial electron and positron. The generator can be used 
for simulation of both tagged and untagged ISR mea- 
surements. For the processes e^e~ — >■ /i^/i^7, e^e^ 
7r+7r~7, and e+e" — >■ K^K"^, NLO FSR radiative cor- 
rections are implemented. In particular, a hard ISR pho- 
ton can be accompanied by emission of a soft photon 
from the final state. 

For all the hadronic states except the two-body 2tt and 
KK as well as 7r"'"7r~7r°, the structure of the electromag- 
netic hadronic current entering the matrix element of the 
process e"'"e^ — > hadrons is model dependent and the ob- 
ject of a study by itself. This model dependence is the 
second source of the theoretical uncertainty. For most of 
the measurements of multihadron cross sections its con- 
tribution significantly exceeds the 0.5-1.0% uncertainty 
of the radiator function. To estimate the model uncer- 
tainty, the distributions of hadrons in data are compared 
to the corresponding simulated distributions. Usually, 
the difference between the detection efficiency obtained 
with different models of the hadronic currents is taken as 
an estimate of the model uncertainty. 



II. EXPERIMENTAL TECHNIQUES 

A. Tagged and untagged ISR 

There are two approaches for studying ISR events. In 
the first approach, the untagged one, detection of the ISR 
photon is not required, but all the final hadrons must be 
detected and fully reconstructed. The ISR events are se- 
lected by the requirement that the recoil mass against the 
hadronic system be close to zero. The mass dependence 
of the detection efficiency for the process e+e^ 7r+7r~7 
at 2£'o = 10.58 GeV is shown in Fig. [51 The efficiency is 
calculated with the Phokhara event generator in leading- 
order mode. The detector acceptance for charged pious 
is assumed to be limited by the condition 30° < 9 < 150° 
which corresponds to the polar angle coverage of the 
BABAR detector. The solid curve in Fig. [5] represents 
the efficiency for the case of untagged ISR photons. For 
two-pion masses below 3 GeV/c^ the detection efficiency 
is about 10% and changes very slowly with mass. At 
these relatively low invariant masses, pious are produced 
in a narrow cone around the vector opposite to the ISR 
photon momentum and therefore can be detected only if 
the ISR photon is emitted at a large angle. The dotted 
curve in Fig. [B] represents the detection efficiency for the 
case of a tagged ISR photon. The photon polar angle is 
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FIG. 6 The detection efficiency for the process e"'"e~ — >■ 
7r^7r~7 at 2Eo = 10.58 GeV as a function of the 2tt invariant 
mass for untagged (solid curve) and tagged (dotted curve) 
ISR photons. 



FIG. 7 The cos Op depen dence of the detection efficiency for 
the process e'^e~ — >■ pp7 (|Aubert et am2006al ). where 6„ is 
the proton angle measured in the pp rest frame with respect 
to the ISR photon direction. The horizontal hne indicates the 
detection efficiency averaged over cos 6p. 



required to be in the range from 30° to 150°. It is seen 
that tagged and untagged efficiencies are very close in 
the mass range below 3 GeV/c^. For higher masses the 
small-angle ISR begins to contribute to the untagged effi- 
ciency leading to its rapid increase, whereas the efficiency 
for the case of a tagged ISR photon varies insignificantly. 

At i3-factories the untagged approach is used for mea- 
surements of exclusive cross sections for masses of pro- 
duced hadronic systems above 3.5 GeV/ c^. The untagged 
detection efficiency is very sensitive to the angular distri- 
butions of the final hadrons. Therefore this approach 
is suitable for the measurement of hadronic processes 
with well defined dynamics, for example, e^e^ — )■ DD or 
e+e" — >■ D*D. For multihadron final states this strong 
sensitivity to hadron angular distributions can lead to a 
sizeable systematic uncertainty of the measurement. 

All measurements of exclusive cross sections of e^e^ 
annihilation into light hadrons at B-factories were per- 
formed using the tagged approach. In contrast to the 
case of untagged ISR, the efficiency for events with a 
detected photon depends weakly on the angular distri- 
butions of the final hadrons. As an example, the angu- 
lar dependenc e of the detect i on effic iency for the process 
e+e^ — > ppj l Aubert et all l2006al) is shown in Fig. [71 
where Op is the proton angle in the pp rest frame mea- 
sured with respect to the ISR photon direction. This 
advantage of the tagged ISR approach allows one to mea- 
sure the cross section for multihadron final states with a 
relatively small model uncertainty. 

Since ISR photons are emitted predominantly along 
the beam axis, in untagged ISR measurements the ad- 
ditional condition that cos6'.y is close to ±1 can be 
used, where 6-y is the polar angle of the momentum 
recoil against the hadronic system in the e~^e~ cm. 



frame. In particular , in Refs. ( Aloisio et all 120051 : 



lAmbrosino et 'all. l2009l) the condition 9^ < 15° or 9^ > 
165° is used to select e+e" — 7r+7r~7 events at the 
0-factory. This condition allows to significantly reduce 
background from the decay (/)—>■ 37r and almost com- 
pletely remove the FSR background, i.e., e~^e~ 7r"'"7r~7 
events with the photon emitted from the final state. It 
should be noted that the FSR contribution related to ra- 
diation by pions is negligible in B-factory experiments 
due to the smallness of the pion electromagnetic form 
factor at s = 112 GeV^. At this energy, the structure- 
dependent contribution, for example, of the processes 
e+e~ /o7 and e+e" — >■ /27 is also expected to be 
small. Theoretical estimations for the cross sections 
of these processes at large s are absent in literature. 
An estimate was made for the process e'*~e~ — )■ pp in 
Ref.(Aubert, 2006a). The FSR contribution (including a 
structure-dependent part) was found to be less than lO^'^ 
for the pp mass below 4.5 GeV. The detection efficiency 



for the process e 



~7 at 2Eo = 1.02 GeV with 



the condition on 9^ described above is shown in Fig. [S] 
The pion polar angles are required to be in the range 50°- 
130°. Due to this restriction the detection efficiency falls 
rapidly with decre asing 27r mass . The untag ged approach 
was u sed in Refs. ( Aloisio et all 120051 : lAmbrosi no et al\ 



|2009( ) to measure the e+e" — )■ tt+tt" cross section in 
the mass range from 0.592 to 0.975 GeV. The tagged 
approach allows one to access the near-threshold mass 
region. The detection efficiency for 7r+7r~7 events with 
a detected photon (50° < 6*^, < 130°) is shown in Fig. [8] 
by the dashed curv e. This selection was also used in the 
KLOE experiment ([Ambrosino et aLl . l201ltlMulleij . l2009l ) 
and allowed to reduce the lower mass boundary for the 
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FIG. 8 The mass dependence of the detection efficiency for 
the process e^e~ — >■ 7r^7r~7 at 2Eo — 1.02 GeV for two 
selections, untagged (6^ < 15° or O-y > 165°) and tagged 
(50° < 9^ < 130°), shown by the solid and dashed curves, 
respectively. The pion polar angles range from 50° to 130°. 



FIG. 9 The mass dependence of the pp mass resolution ob- 
taine d from MC simulatio n for the process e^e~ — > pp-y in 
Ref. (|Aubert et o^.l.l2006a^ . The curve represents the result 
of a polynomial fit. 



cross section measurement from 0.592 to 0.316 GeV. 



B. Hadronic mass resolution and mass scale calibration 

The detector resolution on the hadronic invariant mass 
and the accuracy of the mass scale calibration are impor- 
tant experimental parameters for the ISR cross section 
measurements. 

The mass resolution am is usually determined using 
MC simulation as RMS of the (mmcas — 'm'tmc) distribu- 
tion, where mmcas and mtiuc a-re the measured and gen- 
erated invariant masses, respectively. The experimental 
value of the mass resolution can be extracted from the 
fit of the measured line shape of a narrow resonance, for 
example, J/ip- 

In general, the invariant mass can be represented as a 
sum of the two terms: mmeas = S]^™^ -I- Am(pi,p2, . . .), 
where rui are masses of stable hadrons produced in the 
process under study, and Am is the term depending on 
the final particle momenta pi. The mass resolution am 
is determined by the precision of the measurement of the 
momenta of the charged hadron tracks and photons from 
7r° decays. Since Sim^ has no sizeable spread, and the 
Am term and its uncertainty are minimal near thresh- 
old and grow with the mass increase, it is expected that 
am also increases with mass. As an example, the mass 
resolution versus the p roton-antipro ton mas s for the ISR 
process e+e^ — > pp^ ( Aubert et ali l2006al ) is shown in 
FigEl 

At i?-factories the mass resolution for multihadron sys- 
tems consisting of light quarks varies from 4-7 MeV/c^ 
at the mass of 1.5 GeV/c^ to 6-11 MeV/c^ at 3 GeV/c^; 



the worse values are for hadron states with neutral pi- 
ous. The hadronic cross sections in the mass region be- 
tween the (j)- and J/^-meson resonances do not contain 
structures with a width comparable to the detector reso- 
lution. The 25-MeV/c^ mass bin was chosen for a study 
of most of the processes with light hadrons. With such 
a bin size the distortion of the mass spectrum shape be- 
cause of resolution effects is small. A smaller bin size 
was used for analyses of the processes e^e~ — > ppj and 
e~^e~ 7r+7r^7. For the former, it is important to 
study a near-threshold enhancement in the mass depen- 
dence of the proton electromagnetic form factor. The 
good pp mass resolution for masses below 2 GeV/c^ (see 
Fig|9|) allows to measure the cross s ection in this region 
with the 5 MeV/c^ mass bin C Aubert et a/.l . [20"06a ) . The 
e+e~ — > TT+TT^ cross section near the p-meson peak was 
measured in the BA BAR experiment wit h the mass in- 
terval of 2 MeV/c^ (jAubert et all l2009al) . which is sig- 
nificantly smaller than the tt"*"??" mass resolution (about 
6 MeV/c^ at the p peak). The unfolding of resolution 
effects from the high-statistics (about one-half million 
events) mass spectrum was performed with the proce- 
dure described in Ref. ( Malaescul . l2009l) . The procedure 
uses a mass-transfer matrix that gives the probability 
that an event with true mass in an interval i is recon- 
structed with mnioas in interval j. The transfer matrix is 
usually obtained using MG simulation and corrected to 
take into account a difference in the resolution between 
data and simulation. 

The measure ment of the e"'"e~ — » 7r"*" 7r~ cross section 
at the (/)- factory (jAmbrosino et aLl . [20091) with the KLOE 
detector was performed with the 0.01 GeV^ step in the 



squared mass s' 



corresponding to a mass bin 



width of 6.5 MeV/c near the p peak. The mass reso- 
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lution of the KLOE detector is about 1.3 MeV/c^ at the 
p mass. The resolution effects are substantial only in the 
mass region of the uj-p interference. For comparison with 
theory, these effects were removed by unfolding the mas s 
spectrum using the Bayesian method ( D'Agostini llQQSl) . 

For the J/'ip and ^{2S) produced in ISR processes the 
observed line shapes are fully determined by the detec- 
tor resolution. In this case better mass resolution leads 
to the larger signal-to- background ratio. Fo r the process 
e+e" 2(7r+7r-)7r°7 (jAubert et all l2007cl ) in the mass 
region of the J/ij^ and tp{2S) mesons discussed in Sec- 
tion V, the value of the mass resolution obtained from 
the fit to the J/tp spectrum is about 9 MeV/c^, in good 
agreement with MC simulation. 

For the final states containing charmed and charmo- 
nium mesons ( J/V'7r"'"7r~, DD, ...), the typical resolu- 
tion in the 4-5 GeV/c^ mass range is about 5 MeV/c^. 
The corresponding cross sections were measured with the 
20-25 MeV/c^ mass bin. For these final states the influ- 
ence of the limited mass resolution on the cross section 
measurement is negligible. 

The precision of the absolute mass scale calibration can 
be tested by comparison of the measured mass values for 
known resonances with their nominal values. For many 
multihadron states (see Sec. V) the mass calibration 
is performed at the J/tp m ass. The difference be tween 
the measured and nominal ([Eidelman et al I I2004D J/V^ 
masses is found to be less than 1 MeV /c^ (see, for exam- 
ple, Refs. (lAubert et adl2007cl[2008bl) 1. For the Stt final 
state the mass s cale shift was determ ined at the uj- and 
0- meson masses (jAubert et adl2004bl ): m„ - to^"™'"'*' = 
-(0.2 ± 0.1) MeV/c2 and - m^"™^^ = -(0.6 ± 0.2) 
MeV/c^. We conclude that for the measurements of 
hadronic cross sections at B-factories the mass scale is 
defined with a relative accuracy better than or about 
5 X 10"4. 



C. ISR luminosity 

It is clear that a radiation of a hard photon significantly 
decreases the cross section, so the ISR technique can be 
efficient at high-luminosity colliders only. To compare the 
effectiveness of the ISR method for the measurement of 
hadronic cross sections with direct e~^e~ experiments, it 
is useful to introduce the concept of ISR luminosity. The 
mass spectrum for the ISR process e+e~ — )■ is ex- 
pressed in terms of the ISR differential luminosity dL/dm 
and the Born cross section for the process e+e~ — > X as 

£=e(m)(l + ^(m)).o(™)£, (9) 

where 1 + S{m) = W{m)/WQ{m) is the radiative cor- 
rection factor discussed in Sec. II. CI The ISR luminos- 
ity is proportional to the total integrated luminosity L 
collected in an experiment and the lowest-order radia- 
tor function given by Eq.© or Eq.® depending on the 
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FIG. 10 The mass dependence of the ISR differential lumi- 
nosity multiplied by the detection efficiency. The solid curve 
shows the edL/dm for the B-factory {2Eo = 10.58 GeV, 
L = 500 fb~^, tagged ISR photon), while the dashed curve 
shows the same function for the ^-factory {2Eo = 1.02 GeV, 
L = 240 pb~^, untagged ISR photon). The histogram rep- 
resents integrated luminosities coll ected in direct e"'"e~ ex- 
periments with the SND detector (|Achasov et all |2002| ) at 
the Novosibirs k VEPP-2M collide r (below 1. 4 GeV/c^), and 
with the DM1 l|Bisello et aLlligSlD and DM2 (|Antonelli et all 
11992] ') detectors at the Orsay DCI coUider (above 1.4 GeV/c^). 



angular range used for determination of the detection ef- 
ficiency e(r7i): 



— = Wo TO L. 

dm s 



(10) 



The mass dependence of the ISR differential luminos- 
ity multiplied by the detection efficiency for the BABAR 
experiment is shown in Fig. [10] for masses below 2.2 
GeV/c^. The detection efficiency used was calculated in 
Sec. III. Al for the process e+e" — t- 7r+7r~7 with a tagged 
ISR photon. The integrated luminosity is taken to be 500 
fb^^. The dashed curve in Fig. [TOl shows the same quan- 
tity calculated for the KLOE experiment with the inte- 
grated luminosity of 240 pb~^ and detection efficiency 
taken for the case of an untagged ISR photon (Fig. [5]). 
The luminosity of 240 pb~^ was used in the recent mea- 
surement (Ambrosino et all 12009) of the e+e~ — )■ tt+tt" 
cross section in the 0.592-0.975 GeV/c^ mass range. The 
total integrated luminosity collected by the KLOE is 
about an order of magnitude larger, 2.5 fb~^. The KLOE 
ISR luminosity is shown only up to 0.92 GeV/c^. It in- 
creases sharply and reaches 21 pb~^ at 0.975 GeV/c^. 
It should be noted that the BABAR measuremen t of the 
e+e" TT+vr^ cross section ( Aubert et a^.l . [2009al) is also 
based on a part of the recorded data corresponding to 232 
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FIG. 11 The mass dependence of the ISR differential lumi- 
nosity multiplied by the detection efficiency for experiments 
at the B-factory {2Eo = fO.58 GeV, L = 500 fb~\ untagged 
ISR photon) in the charm production mass region. 



fb""'^. The histogram in Fig. [TOl shows the distribution of 
the integrated luminosities cohected in some direct e^e^ 
experiments. At masse s below 1.4 GeV/c^ t he statistics 
of the SND experiment ( Achasov et a/.l . [2002t ) recorded at 
the VEPP-2M collider is presented. This is the largest 
integrated luminosity collected in this mass region in a 
single experiment. The mass bin 1.0-1.1 GeV/c^ does not 
include about 13 pb^^ taken by SND in vicinity of the 
0-meson resonance. The significant part of the statistics 
from the 0.7-0.8 GeV/c^ mass interval is collected in the 
w-meson mass window 0.76-0.80 GeV/c^. In the cm. 
energy range 1.4-2.2 GeV/c^ the experiments with the 
largest statistics are DM1 and DM2 at the Orsay DCI 
e+e" collider. The histogram at to > 1.4 GeV/c^ shows 
a sum of the integrated luminosities collected with these 
detectors. 

At low masses of the hadronic system the data sam- 
ples of ISR events currently available at i?-factories ex- 
ceed the statistics collected in conventional e+e" experi- 
ments, especially at masses below 0.7 GeV/c^ and above 
1.4 GeV/c^. The ISR luminosity of the </)- factory in- 
creases very rapidly with mass. For masses below 0.8 
Ge V/c^ the lumin osity currently used for ISR analy- 
sis ([Ambrosino et al . , 2009) is comparable with that col- 



lected in direct e+e" experiments. For higher masses it 
exceeds both BABAR and e+e" luminosities. 

The ISR luminosity for the mass region of charm 
production is presented in Fig. [TT] It corresponds to 
the 500 fb~^ integrated luminosity collected at 2Eo — 
10.58 GeV/c^ and is multiplied by the detection effi- 
ciency calculated for the case of an untagged ISR photon 
(Fig. [B]). The ISR luminosity in this mass region sig- 
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FIG. 12 The de tection efficiency for the process e'^e~ 
ppy (jAubert et a l.. 2006a) as a function of the pp invariant 
mass. The curve represents the result of a polynomial fit. 



nificantly exceeds the integrated luminosity collected in 
direct e^e~ experiments including t he re cent CLEO-c 



energy scan ( Cronin-Hennessv et al\ . [20091 ) . 60 pb~^ at 
twelve points between 3.97 and 4.26 GeV. 

Thus, the current data samples of ISR events produced 
at the B- and ^-factories are larger than those produced 
directly in e+e" collisions for all masses of interest ex- 
cluding the regions near the narrow resonances (w, 0, 
J/f/;, il]{2S)). For masses above 1.4 GeV/c^ this allows to 
significantly improve accuracy of the measurements of ex- 
clusive hadronic cross sections. In the mass region below 
1.4 GeV/c^ the results obtained with the ISR method are 
comparable to rather precise direct e+e" measurements. 



D. Comparison with e+e scan 



The ISR technique offers some advantages over con- 



ventional 



measurements. One of them is that the 



entire hadronic mass range is accessible in one experi- 
ment. This allows one to avoid relative normalization 
uncertainties which inevitably arise when data from dif- 
ferent experiments, or from different machine settings in 
one experiment, are combined. 

The ISR measurements with a tagged photon have ad- 
ditional advantages. In many cases, particularly for final 
states with low invariant mass of the produced particles, 
the hadronic system is coUimated along the direction op- 
posite to the ISR photon. Therefore, the detection effi- 
ciency has low sensitivity to hadron angular distributions 
in the hadronic-system rest frame. In Fig|7] the angular 
dependence of the dete ction efficiency is sho wn for the 
process e"*"e~ — >■ pp7 (jAubert et all l2006at) . The an- 
gular dependence is close to uniform. This reduces the 
model dependence of the cross section measurement due 
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to the unknown relation between the values of the proton 
electric and magnetic form factors, and significantly fa- 
cilitates data analysis. Note that in conventional exper- 
iments at e^e~ or pp colliders the detector acceptance 
for the final pp or e'^e~ systems falls to zero when cos9p 
approaches ±1. 

For ISR events the final hadrons have non-zero mo- 
menta at the production threshold and are there- 
fore detected with full efhciency. In FigH^] the de- 
tectio n efficiency f o r the process e+e" — ?> ppj pro- 



cess ( Aubert et all . l2006al) is shown as a function of 
the pp invariant mass. No strong variation of the ef- 
ficiency with mass is observed, while in direct e^e~ 
measurements the detection efhciency vanishes at the 
threshold because of the low momenta of the pro- 
duced particles. This feature of ISR hadron produc- 
tion was successfully used at BABAR for t he mea sure- 
ments of the e" ^e~ — >■ pp 



(lAubert et all l2006al ) and 



2009a|) 



cross sections m 



( Aubert et all 
near-threshold mass regions. 

For the measurement of the e+e~ — >■ tt+tt" cross sec- 
tion, particle identification plays a crucial role. In the 
ISR process e+e^ — > 7r+7r^7 at i?-factories most of the 
final pions have momenta larger than 1 GeV/c. For such 
pion momenta a good 7r//i/e separation is provided which 
allows one to almost completely remove the e~^e~ — > ■ 
e"'"e~7 and e+e~ — >■ /i+/x~7 backgrounds ( Aubert et all 
l2Q09aD . This is in contrast with direct e~^e measure- 
ments ([Achasov et all 120061: lAkhmetshin et all l2004al 



I2007D in which it is difficult to separate e"* 
and e'^e^ — > events in the most interesting p- 

meson mass region (0.60-0.95 GeV). As a result, a sum 
of the cross sections is measured. The contribution of 
the process e^e~ — fJ.~^fJ-~ is then subtracted using its 
theoretically calculated cross section. This leads to an 
increase of statistical and systematic errors of the mea- 
surement. 

It should be noted that the advantages of tagged ISR 
discussed above (weak mass and angular dependences of 
the detection efficiency) are completely absent for un- 
tagged ISR. In this case the mass and angular depen- 
dences are even stronger than those for events of direct 
e'^e^ annihilation. 

A disadvantage of ISR is that the mass resolution and 
absolute mass scale calibration are much poorer than the 
beam energy spread and the accuracy of the beam energy 
setting in direct e^e~ annihilation experiments. The in- 
fluence of the resolution effects on the ISR measurement 
is discussed in Sec. III.Bl 

The main disadvantage of the ISR measurements is 
presence of a wide spectrum of background processes dif- 
ferent from those in direct e^e~ experiments. For ex- 
ample, in e+e~ annihilation the main background pro- 
cess for e+e^ 
a lost 



For the ISR process e+e 



with the Stt mass in the range msTt ± Ato/2, this back- 
ground corresponds to the contribution of the process 
e~^e~ — 7r^7r~7r''7r''7 with the An mass in the same range 



msTT ± Am/2. The presence in ISR of Attj events with 
arbitrary masses, which may, in particular, be out of the 
msTT ± Am/2 range, greatly increases background. 

At the (/>-factory and in future ISR measurements at 
the tau-charm factory in Beijing the background from 
FSR processes should be taken into account when the 
ISR photon is detected. The FSR contribution for the 
e+e" — ?> 7r+7r~ measurements at KLOE is calculated 
with the PHOKHARA generator, which models FSR for 
pions using scalar QED, and also takes into account the 
radiative 4> decays to tt+tt''' via the /o (980)7 ^-nd pir in- 
termediate states. The pion electromagnetic form fac- 
tor used in the generator is obtained from a fit to the 
e+e^ — > TT+TT^ experimental data. In the case of the 
tau-charm factory, experimental information on exclu- 
sive hadronic cross sections in the energy range from 3.0 
to 4.5 GeV obtained at B-factories can be used to esti- 
mate the FSR contribution. Additional theoretical input 
is required to estimate structure-dependent FSR. 

Another background source is the non-ISR process of 
e+e^ annihilation into hadrons containing a high-energy 
7r°. In particular, the events of the process e~^e~ Xtt° 
with an undetected soft photon or merged photons from 
the 7r° decay may almost completely imitate the e'^e^ -^• 
events. This background is usually subtracted sta- 
tistically using for normalization selected e+e~ — )• Xtt° 
events with a reconstructed 7r°. In tagged ISR mea- 
surements at i?-factories the process e+e^ — ?> Xtt" be- 
comes the dominant background source at relatively high 
masses, about 2 GeV/c^. It limits the mass region for ISR 
studies of light hadrons to masses below 4.0-4.5 GeV/c^. 

In ISR measurements with an untagged ISR photon, 
the background from e+e^ — >■ Xtt^ can be significantly 
suppressed by requiring that the missing momentum in 
an event be directed along the beam axis. For untagged 
ISR, the main sources of background are ISR processes 
and two-photon processes e"*'e~ — >■ e+e~7*7* — ^ e~^e~ X 
in which initial electron and positron are scattered pre- 
dominantly at small angles. The latter background can 
be suppressed by a condition on the missing mass, which 
should be close to zero for ISR events and has a wide 
distribution for two-photon events. 

Background suppression and subtraction are the main 
sources of the systematic uncertainty on ISR measure- 
ments. 



E. Colliders and detectors using ISR 



ISR processes were studied in many e^e~ experiments 
either as a source of useful physical information or as 
a source of background. For example, possibly the first 



study of the process e^e 



"TT 7 was performed more 



than 20 years ago with the ND detector at the VEPP-2M 
coUider (jPohnskv et all Il99ll: IVasserman et all Il988[) . 



In this work, the FSR process e+e" — )■ /: 
measured with the ISR process e~^e~ p"f ^ tt^'tt^j 
studied as a main source of background. Many inter- 



TT^TT 7 was 
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esting ISR studies have been performed with the CLEO 
detector, see, e.g., Ref. ()Adams et "oZI. I2006D . Below we 
give a brief description of only three detectors: BABAR, 
Belle, and KLOE, which made a great contribution both 
to development of the ISR technique and ISR measure- 
ments of hadronic cross sections. 



F. PEP-II and BABAR 

The PEP-II S-factory at SLAC is a two-ring 
asymmetric-energy e^e^ collider with energies of 9 GeV 
for the electron and 3.1 GeV for the positron beam, oper- 
ating at the cm. energy of 10.58 GeV , at th e maximum 
of the T(4S') resonance ( Seeman et The max- 

imum luminosity achieved at PEP-II was slightly over 
1034 cm-2 s-\ The principal goal of the PEP-II B- 
factory and the BABAR detector is studies of CP viola- 
tion in the _B-meson system. 

The BABA R detector (F i giT^ is described in de- 
tail elsewhere ( Aubert et all [20021) . Final states with 



charged particles are reconstructed in the BABAR track- 
ing system, which comprises a five-layer silicon vertex 
tracker (SVT) and a 40-layer drift chamber (DCH) oper- 
ating in a 1.5-T axial magnetic field. The vertex position 
is measured by the SVT with the accuracy of 50 /xm. 
The momentum resolution for 1 GeV/c charged tracks 
is cTpt/pt = 0.5%. Charged-particle identification is pro- 
vided by an internally reflecting ring-imaging Cherenkov 
detector (DIRC), and by energy loss measurements in 
the SVT and DCH. The hard ISR photon and photons 
from 7r° decays are detected in a CsI(Tl) electromagnetic 
calorimeter (EMC). The energy resolution for 1 GeV pho- 
tons is about 3%; the angles of photons are measured 
with the 4 mrad accuracy. Muons are identified in the in- 
strumented fiux return (IFR) of the solenoid, which con- 
sists of iron plates interleaved with resistive plate cham- 
bers. 

Experiments at the PEP-II colhder with the BABAR 
detector were carried out from 1999 to 2008. The to- 
tal integrated luminosity is close to 530 fb^^. The ISR 
studies at BABAR started in 2001. The ISR research 
program includes a study of the light hadron production 
with a tagged ISR photon and charm and charmonium 
studies with an untagged photon. 



G. KEKB and BELLE 

The KEK B-Factory, KEK-B, is an asymmetric-energy 
(similar to PEP-II) e+e" collider with the 8-GeV electron 
and 3. 5- GeV positron beams and th e maximum lu minos- 
ity of 2.1 • 10^4 cm-2 g-i (Kuroka wa et all l2003l ). The 
main physical goal of this project is to perform a detailed 
study of B-meson prop erties, in particular, C P-violation. 

The Belle detector ()Abashian et all |2002|) (FigH!]) is 
configured inside a 1.5 T superconducting solenoid. The 
_B-meson vertices are measured in a three-layer dou- 




FIG. 13 View of the BABAR detector. 
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FIG. 14 Side view of the Belle detector. 



blesided silicon vertex detector with about 50 fim im- 
pact parameter resolution for 1 GeV/c momentum track 
at 6 ~ 7r/2. Track momenta are measured in a 50-layer 
wire drift chamber with a 0.4% momentum resolution at 
1 GeV/c. Particle identification is provided by dE/dx 
measurements in the drift chamber, aerogel Cherenkov 
counters, and time-of-flight counters placed outside the 
drift chamber. Electromagnetic showers are detected in a 
CsI(Tl) calorimeter located inside the solenoid coil. The 
energy resolution is 2% for 1-GeV photons. An iron flux- 
return located outside the coil is instrumented to detect 
KL-mesons and identify muons. 

Experiments with Belle started in 2000 and stopped 
in 2010. The Belle integrated luminosity reaches 1000 
fb^^. The ISR experiments are mainly devoted to the 
production of charm and charmonium hadronic states 
with mass above 4 GeV/c^. ISR analysis of light mesons 
is in progress. 
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FIG. 15 KLOE detector (jMulleJ. |2009| ). The polar-angle 
regions used to select tagged (50° < 9-^ < 130°) and untagged 
(61^ < 15° or 6^ > 165 ) ISR, events are shown. 



H. DAcE-NE and KLOE 



DAi>NE, 



the 



Frascati 



-factory 



( Franzini and Moulsonl . l2006[ ). is in operation since 
1999. The main goal of the DA$NE project is a study 
of neutral and charged kaons, intensively produced at 
the energy corresponding to the maximum of 0(1020) 
resonance. Similar to PEP-II and KEK-B, DA$NE uses 
two separate rings for storing electron and positrons, 
but beams have equal energies. The DA'l'NE design 
luminosity is 5 • 10^^ cm~^ s~^. 

KLOE (jFranzini and Moulsonl . |2006( ) (Fig lT5|) is the 
main DA$NE detector. The detector consists of a large- 
volume drift chamber (DC) surrounded by a hermetic 
electromagnetic calorimeter (EMC). A superconducting 
coil provides an axial magnetic field of 0.52 T. In or- 
der to reduce neutral kaon regeneration and charged- 
particle multiple scattering, the gas mixture of 90% he- 
lium and 10% isobutane is used in the DC. Charged- 
track momenta are measured with the <Jplp = 0.4% 
accuracy. The lead-scintillation fiber calorimeter pro- 
vides the energy resolution for electromagnetic showers 
of (7e/E = 5.7%/-\/i?(GeV), and the time resolution of 
at = 54ps/^£;(GeV)014Ops. 

The total integrated luminosity accumulated with 
KLOE is about 3 fb~^. The only, but very important, 
ISR process studied at KLOE is e+e" 
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III. PRODUCTION OF LIGHT QUARK MESONS 
A. Overview 

As already mentioned in the Introduction, e+e^ anni- 
hilation into hadrons at cm. energies below 2 GeV plays 
a very important role in many fundamental problems of 
particle physics. In particular, knowledge of its total 
cross section is mandatory for the calculation of the muon 
anomalous magnetic moment in Standard Model. For 
many years only e^e^ scan experiments provided infor- 
mation on this reaction and determined t he uncertainty 
of the SM prediction of the muon anomaly ([Pavier et all . 
I2003al lbh. Main information on light vector mesons has 
been also obtained in such measurements. Unfortunately, 
the collected data samples were not sufficient for a precise 
determination of parameters of excited vector mesons. 
Recently, due to a very high luminosity of the e^e~ fac- 
tories, DAFNE, KEK-B, and PEP-II, the ISR technique 
became a powerful tool for an independent study of e^e^ 
annihilation at low energies. 

The KLOE collaboration used the ISR method at 
the (f) meson energy to study the reaction e+e^ — > 
TT+TT " and measu r e the pion electromagne t ic form fac- 
tor (lAloisio et all 120051 : lAmbrosino et all l2009l l201ll: 
lMulleil . l2009ir Recently, results on this pr ocess were also 
report ed by the BABAR collaboration ( Aubert et all . 
l2009aD . 

of high-multiplicity final states were 
BABAR: tt+tt-tt' 



A variety 
studie d at 
l2n04bD . 

2{K+K-) 



2(7r+7r~7r° 



!(7r+7r-). 



l2006bD . 

K+K 
and K+K 



(jAubert et al 
and K+K- 
2(7r+7r-)7rO, 2 



TT+TT tt" (jAubert et all . 
Tr+TT~ K + K~ and 
ioolE), 3(7r+7r-) 
'2(7r+7r~) (lAubert et all . 



-1] 



TT+TT )ri, K+K TT+TT tt'^ and 



2007cD 
' 2007bl). 



K+K- 
K+K- 



Aubert et all 
-"tt" (Aubert et all 
and K+K-K+K- 
K+K% Tr+, K+K- Tr° and K+K'r] (jAubert et all . 
l2008bh . The fi nal K+K' 'TT+TT state was also investi- 
gated by Belle (jShen et aLl . l2009l ). 



(lAubert et ail. l2007bl) 



Studies of the exclusive channels of e+e~ annihilation 
listed above allow to determine such fundamental param- 
eters as mass, width and leptonic width of various vector 
mesons. In addition to the low-lying resonances, such as 
the p, w and (j), where ISR studies can independently pro- 
vide meaningful and competitive information, they are 
indispensable for a much more precise than before inves- 
tigation of the excited vector states. 

Moreover, detailed analysis of the dynamics shows that 
in many cases a multiparticle final state can be reached 
via different intermediate mechanisms. For example, four 
pious can be produced via W7r°, ai(1260)7r, p^fo, .... In 
the following sections we show a complexity of the in- 
ternal substructures observed in some channels, which 
are often used to extract parameters of the resonances 
involved in the substructures. 

In general, amplitudes corresponding to different inter- 
mediate mechanisms interfere affecting the energy and 



50 
45 
40 
35 
30 
25 
20 
15 
10 
5 




> KLOE09 
ACMD 



0.1 0.15 0.2 0.25 0.3 0.35 



\ A CMD 
★ SND 




KLOE09 



I .... I 



(M;?^' [GeV'] 



I .... I 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



0.1 
0.05 


-0.05 
-0.1 




(|fJc,s-IFJk)/|fJk 



4ri- 



*-4 



* CMD-2 

* SND 



0.4 0.5 0.6 0.7 0.8 0.9 



1 



FIG. 16 Top: The pion form factor obtained by KLO E in the 
react ion e+e^ — >■ TT+TT~y with a tagged ISR photon (| Mulled . 
l2009l ). Bottom: Relative difference between the KLOE re- 
sult with an untagged ISR photon (jAmbrosino et all. |2009|) 
and t he direct e+ e~ measurements by SND (jAchasov et all 
|2006| ) and CMD-2 (jAkhmetshin et allhoO'^ . The dark (light) 
band indicates the KLOE uncertainty (statistical and sys- 
tematic errors combined in quadrature). For the SND and 
CMD-2 data, the combined statistical and systematic errors 
are shown. 



angular distributions of the final particles. This inter- 
ference should be taken into account to avoid additional 
systematic errors. 

Unless otherwise stated, all cross sections in the fol- 
lowing sections are corrected for effects of initial-state 
radiation only. Neither final-state radiation nor vacuum 
polarization corrections have been applied. 



B. 



e e — >■ TT TT 



The reaction e+e~ — > tt+tt~ was relatively well studied 
for cm. energies up to 1.4 GeV in direct e+e~ experi- 
ments. The most precise measurements were perfo rmed 
with t he CMD-2 (jAkhmetshin et all j2004al j2007D and 
SND (jAchasov et all , j2006^ detectors at the VEPP-2M 



14 



0.4 
0.3 
0.2 
0.1 


-0.1 
-0.2 
-0.3 
-0.4 



■ 



▲ Data KLOE09 (Vs = 1 GeV) i+i»*X 
T Monte Carlo 



I ' ' ' ' I I I ' ' I ' I ' ' I I I I ' I ' I I I I ' ' ' ' I ' ' I I I ' ' ' ' I I I ' 



itive with them in a s ystematic uncertainty. Both 
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FIG. 17 Forward-backward asymmetry in tlie reaction 
e^e~ — > 7r^7r^(7) measured with the KLOE detector (|Mulleij . 
120091 '). 
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FIG. 18 Top: The QED test by the ratio of the e+e" -5> 
fj,'^ fj.~ cross section in data to the theoretical one. Bottom: 
The e'^e~ — >■ ■K'^n~ cross sectio ns measured with the BABAR 
detector ijAubert et d.Ll2009al '). 



collider. The CMD-2 measurements have a systematic 
uncertainty in the 1% range. 

The dominant contribution to this process comes 
from the p(770) meson production. A measure- 
ment of the e'^e~ — )■ 7r"^7r~ cross section in the 
p(770) mass regi on was perf o rmed by KLOE using 
the I S R me t hod (lAloisio et all . l2005l : lAmbrosino et all 
I2009L 1201 ll: iMulled . I2009D . For the first time it 
was demonstrated that the cross section determined 
by this method could have smaller statistical errors 
than direct e~^e~ measurements and could be compet- 



untagged ( Aloisio et al . 20051 Ambrosino et aL ^ 20091 ) 



and tagged (|Ambrosino et ail . 1201 ll : iMulleil . 120091 ) ISR 

7r+7r~7 events were studied with consistent results. 
While the tagged measurement has worse statistical er- 
rors and an additional source of the systematic uncer- 
tainty due to the FSR contribution, it covers the re- 
gion of small invariant masses inaccessible for the un- 
tagged measurement^The^ measure- 
ment ()Ambrosino et all l201ll: iMuUeil . l2009l) represented 
as a pion electromagnetic form factor squared is shown 
in Fig. fW top) together with the r esults of the direct 
e'^e~ measureme nts with the CMD-2 dAkhmetshin et all 
l2007t ) and SND (|Achasov et all [20061) detectors. Com- 
pariso n of the more pr ecise untagged KLOE measure- 
ment (| Ambrosino et aLl . 12009 ) with the CMD-2 and SND 
data is given in Fig. ITBT bottom) . At invariant masses cor- 
responding to the maximum of the p resonance and its 
high-mass tail the points from direct e~'"e~measurements 
lie systematically higher than those from KLOE. In 
this mass region the difference between the CMD-2 and 
KLOE measurements is definitely larger than their com- 
bined systematic uncertainty. The KLOE systematic 
error includes the experimental (0.6%) and theoretical 
(0.6%) uncertainties. Two main sources of the former 
are tracking and luminosity measurement. The latter is 
determined mostly by the accuracy of the radiator func- 
tion calculated with the PHOKHARA event generator. 
Note that the KLOE Collaboration performed a dedi- 
cated study to validate a calculation of FSR effects us- 
ing forward-backward asymmetry arising from the inter- 
ference between the ISR and FSR amplitudes (Miiller, 
2009). The study showed that the assumption of point- 
like pions works reasonably well and can be used for the 
FSR calculation, see Fig. [T71 

A structure seen at the top of the p- meson resonance is 
due to its interference with the much more narrow a;(782) 
resonance also decaying to tt+tt". Because a;(782) mass 
is known precisely, the position of this structure can be 
used to test the accuracy of the mass scale calibration. 
Unfortunately, neither KLOE nor BABAR (see below) 
report the result of such a test. 

The PEP-11 B-factory also provided a large sample of 
the e+e" -> 7r+7r~7 events (abo ut 530 thousand) and 
the e+e^ — > tt+tt^ cross section ([Aubert et al I l2009aD 
was measured for the e"'"e^ cm. energies up to 3.0 GeV. 
In this experiment the e+e^ — >■ 7r"'"7r~ cross section is 
obtained from the ratio of the tt+tt^ and /i+/i^ mass 
spectra. Due to the normalization to the cross section of 
the theoretically well known process e~'"e~ — > fi~j, the 
measurement becomes much less sensitive to the exper- 
imental uncertainties and to the theoretical uncertainty 
of the radiator function. A comparison of the measured 
mass spectrum for the reaction e+e" — > /i+/i^7 
with the QED prediction is shown in Fig. fTSFtop). The 
data and the prediction are consistent within the esti- 
mated systematic uncertainty of 1.1%, dominated by the 
accuracy of the integrated luminosity measurement. Us- 
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FIG. 19 The e+e~ — >■ 7r"'"7r~ cross se ction above 1 
GeV measured with the BAB AR detector (Wangl 120091 ). 
Compari son with t he CMD-2 l|Aulchenko et all [2005h and 
DM2 (jBisello et all il989. ') measurements is shown. 
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2 dAkhmetshin et all l2004al . |2007| ) and BABAR l|Wang I 
measurements. The band corresponds to the BABAR 
statistical and systematic uncertainties combined in quadra- 
ture. 



The claimed sub-percent level of systematic uncertain- 
ties on the e'^e~ — >■ 7r"'"7r~ measurements can be verified 
by comparison of the results from these very different ex- 
periments. Above we found that the difference between 
the KLOE and CMD-2 measurements is larger than their 
combined systematic uncertainty. Figure shows a rel- 
ative difference between the KLOE and BABAR mea- 
surements. Again the deviations larger than declared 
systematic errors are seen indicating a presence of unac- 
counted systematic uncertainties in one or both exper- 
iments. Comparison between the CMD-2 and BABAR 
shown in Fig. [21] also reveals some non-statistical up 
to 5% deviations both below and above the p-resonance 
maximum. In the whole energy ra nge BABAR d a ta are 
in fair agreement with the SND ( Achasov et all . l2006l ) 
results within experimental uncertainties. 

In Section IVIII we discuss the impact of these measure- 
ments on the problem of the muon anomaly. 



ing the bin-by-bin ratio to the cross section of the pro- 
cess e^e^ — > //^/i^7 one minimizes theoretical uncer- 
tainties and reduces a systematic error at the p peak 
to 0.5% dominated by pion identification and ISR lumi- 
nosity. Previously such a test was performed in e+e" 
scan experiments at the OLYA detecto r in the cm. en- 
ergy range from 640 MeV to 1400 MeV (jKurdadze et all . 
Il984l ) an d at the CMD-2 detec tor from 370 MeV to 
520 MeV (jAulchenko et a/.l . l2006l ) with the achieved pre- 
cision of comparison of 3% and 1%, respectively. 

The measured cross section is shown in 
Fig. [TSl' bottom). For the first time a relatively 
high-statistics measurement is performed for cm. en- 
ergies above 1 GeV. The cross section in this energy 
range shown in Fig. [12] demonstrates some statistically 
significant structures which can be possibly explained by 
the interference between the wide p-like excited states. 
Note that the cross section shown in Fig. [TSjis bare and 
includes FSR effects. 



C. 



e e — > TT TT TT 



A study of the three-pion production i n the ISR pro- 
cess w as reported by BABAR in Ref. ( Aubert et all . 
l2004bD . The three-pion mass distribution for the e+e — >■ 
7r+7r~7r°7 reaction shown in Fig. [22] is dominated by the 
well known a;(782), (/)(1020), and J/ijj resonances. For 
the cij(782) and ^(1020) resonances they determine the 
product of the leptonic width and the branching frac- 
tion to three pions consistent with other measurements 
and having comparable accuracy. Large data samples 
make possible the observation of two structures in the 
1-2 GeV/c^ mass region (see Fig. The cross sec- 

tion below 1.4 GeV is in agreem ent with the SND mea- 
surement (jAchasov et all |2002| ). but at higher energies 



a larg e deviation from the DM2 results ( Antonelli et all . 
Il992f ) is observed. The cross section in this region is 



fitted (see inset in Fig. [^5)) assuming the p resence of two 
excited cj-like states, cli(1420) and w(1650) (jAmsler et all . 
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FIG. 22 The m(7r+7r"7r°) distribution for the e+e" -5> 
7r"'"7r ~7r°7 reaction mea sured with the BABAR detec- 
tor (|Aubert et all\2QQA\i ). 
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FIG. 23 The e^e~ — >■ 7r"'"7r~7r° cross sect ion measured with 
the BABAR detector (lAubert et a/.Ll2004bl') in the 1-3 G eV/c^ 
range compared wit h the SND |Achasov all l2002f ) (open 
circles) and DM2 l|Antonelh aLl 1 19921 ) ftriangles) data. 
The inset shows the mass distribution fitted with two res- 



I2OO8I ). The parameters of these states are stiU not weU 
determined. In this case they strongly depend on rel- 
ative phases between the corresponding amplitudes and 
their phase differences with the a;(782) and 0(1020) am- 
plitudes. The latter resonances have a much larger decay 
rate to the Stt mode. The obtained parameters of the 
w(1420) and a;(1650) are summarized in Table HIl 

The three-body final state is a relatively simple process 
for a study of hadron dynamics. Its Dalitz plot analysis 
shows that the p{770)'^n^ and p{770)°n^ intermediate 
states dominate at all energies. There is also a small 
contribution of the wtt intermediate state with ui decay 
to TT+vr^. 



D. 
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FIG. 24 The e+e" TT^ cross section measured 

by BABAR (|Aubert et oU (2008bl ) (top). Comparison of 
the BABA R measureme nt with the results of the prev ious 
DM1 ( Buon all Il982l ) and DM2 (|Bisello efail ll991bD ex- 
periments (bottom). 



(top) and comparison of the B ABAR resul t s with 
DM1 (|Buon et all\l98'^ and DM2 (jBisello et aLl . ll991bf ) 
measurements below 2.4 GeV where the previous data 
are available (bottom). The BABAR data are about 10 
times more precise. The "spike" at 3.1 GeV is due to J/ip 
decays to these final states and will be discussed later. 

In the K+K-2j final state the (/)(1020)77 and 0(lO2O)7r 
intermediate states were also observed. The measured 
cross sections for these states, which were not previously 
studied, are shown in Figs. [51] and [571 The e+e^ — >■ 
0(1020)77 is the best channel for a study of the excited 
6 state. The contributions of the w-like states to this 
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ment (|Bisello et all ligQlbl ) (bottom). 
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FIG. 27 The (j!)7r° cros s sectio ns measured with the BABAR 
detector l|Aubert et aUbOOSbl ). 



channel should be suppressed by the OZI rule. 

The reaction e^e~ — >■ (/)(1020)7r° is promising for a 
search for exotic isovector resonances. For ordinary 
isovector states, the (f)Tr'^ decay should be suppressed by 
the OZI rule. The authors perform two fits of the cross 
section. In the first one they assume a single resonance 
and obtain for it mass and width of 1593±32 MeV/c^ and 
203 ± 97 MeV, respectively. These parameters a. r e cora - 
patible with those of the p(1700) (|Amsler et all . |2008[) . 
A somewhat better quality of the fit is achieved if two 
resonances are assumed. The obtained parameters of the 
first resonance are 1570±36±62 MeV/c^ for the mass and 
144±75±43 MeV for the width, i.e., consistent with those 



of th e C(1480) state observed in Ref. ( Bitvukov et al\ . 
I1987D . The mass and width for the second resonance 
are 1909 ± 17 ± 25 MeV/c^ and 48 ± 17 ± 2 MeV, re- 
spectively, compatible with the dip already observed 
in other experime nts, pre dominantl y in mu ltipion fi- 
nal states (.A nt onelli et al\. |l996: Aubert et al.. . .2006bl : 
iBaldini et all . ir988HFrabetti et a/.l . ]2000l) . With the lim- 
ited statistics available at the moment they cannot draw 
a firm conclusion: an OZI-violating decay of the p(1700) 
cannot be excluded. 

Figure [28] shows the Dalitz plots for the K~^K~7t'^ and 
K'^K^n^ final states. It is seen that the KK* {892) and 
i4ri^2(1430) intermediate states give the main contribu- 
tion to the KKtt production. For the KgK^ir^ final 
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FIG. 28 The Dalitz plot distribution for the K+K~tv° (a) and ifgA'=^7r'F final state (b) from Ref. (jAubert et aZ.I . l2008bD . A 
sum over all accessible cm. energies of the hadronic final states is given. 
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TABLE I Parameters of the isoscalar and iso vector res- 
onances obtained in Ref. l|Aubert et alV [2008b) from the 
global fit to the isoscalar and isovector amplitudes using the 
e+e" ^ A'±A'*(892)T,A'0A"±7rT and e+e" ^ (j)Tj cross sec- 
tions. 



state both the neutral K^K*^ and charged K^K*^ com- 
binations are involved. Since the K°K*^ and K^K*^ 
amplitudes are the sum and the difference of the isovec- 
tor and isoscalar amplitudes, respectively, the Dalitz 
plot population for the K'gK^-K^ mode is asymmet- 
ric and strongly depends on isospin composition. From 
the Dalitz plot analysis the moduli and relative phase 
of the isoscalar and isovector amplitudes both for the 
KK*{m2) and KK^(IA3,Q) intermediate states were de- 
termined. The obtained isoscalar and isovector e^e~ — >■ 
KK*{%Q2) cross sections are shown in Fig. [29] (a, b). 

The global fit to the e+e' 0(1020)77 and 

e+e" — >■ K+K'TT^ cross sections, isovector and isoscalar 
if j^* (892) amplitudes, and their relative phase was per- 
formed to determine parameters of the (j) and p excita- 
tions decaying into these final states. The fit results are 
shown in Fig. [29] and summarized in Table [l] The ob- 
tained mass and width of the (j)' and p' are in reason- 
able agreement with the parameters of the p(1450) and 
(^(16 80) resonances meas ured in other experiments (see 
Ref. (jAmsler et al. 1. 12008[ ) for references). The parame- 
ters of the (/)", which is seen in the 0?/ final state, are close 
to those for the F(2175) state observed in the (/)/o(980) 
final state. This state will be discussed in Sec. IIII.FI 



E. 



■K + TX TT + TX , TT^TT 2^" 



The reactions e+e~ — >■ 7r"'"7r~7r"'"7r~ , 7r"'"7r~27r'' have 
the largest cross sections in the energy region above the 
0-meson resonance. Th e y were studied with the B ABAR 
detector (jAubert et all l2005bl: iDruzhininl . l200iD in the 
energy region below 4.5 GeV. Figures [30] and [31] show 
comparison of the BABAR results wit h the previous 
direct e+e" measurements , see R efs. (lAchasov et 



20031: lAich metshin all. l2004bl: iBacci aLl. ll98C ; 



Bisello et ali Il991at ICordier et al\. Il982at [Cosrne et al. 



19791 : iDolinskv et all Il99ll: iKurdadze et all Il988l ) 
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R with 7 = 


0' <P" 
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rfeSjJx* (892) (eV) 
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m_B(MeV) 
rfl(MeV) 
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R with 1 = 1 
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m_R(MeV) 
rH(MeV) 


135 ± 12 
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30 Comparison of the BA BAR results on the 
e+e' — >■ tv+tt' -k+tt' cross section (lAubert et all l2005bD 
with the previous dir ect e^ e~ measu rements (Achasov et al , 
20031: Akh mctshi n et al\. [2004bl: jBacci_ei al., 198C 



Bisello et ali limid ICordier et all Il982al: ICosme et al. 



19791 : IDolinskv et aZ.t[l9'9ll : lKurdadze et aLlll988l ) 



for 



TT TT TT TT 



and R efs. (lAchasov et all 2005 ; 



Akhmetshin et al. ^ 1999; B acci et all 1981 1 Bisello et al. 



Kurdadze et all\l98($) for tt+tt 27r°, in the energy range 



1991al: ICosme et all Il979[ IDolinskv et"al[ Il991 



covered by these measurements. The large difference be- 
tween the data sets from different experiments indicates 
that some previous measurements had large, up to 50%, 
unaccounted systematic errors. The BABAR systematic 
uncertainty on the e+e' — > tt+tt'tt+tt' cross section 
is estimated to be about 5% in the 1-3 GeV energy 
range. For this channel the BAB AR data are i n good 
agreement wi th the recent SND (lAchasov et al. 1, I2003D 
and CMD-2 ([Akhmetshin et all l2004bD measurements 
at the energies below 1.4 GeV. The DM2 and BABAR 
data are also in reasonable agreement. 

For the 7r+7r^27r° channel the BABAR results are 
still preliminary. The estimated systematic uncertainty 
changes from 8% in the maximum of the cross section to 
10% at 1 and 3 GeV. 

At energies below 1.4 GeV the BABAR cross 
sections agree well with t he resul ts of the re- 
cent S ND (lAchasov et al. 

L I2003D and older 
OLYA ( Kurdadze et al!. '1986*) measurements, but 
not with the ND ( Dolinsky et al., 1991) and CMD- 
2 ( Akhmetshin et ali Il999f ) cross sections that may 



be affected by large unaccounted systematic errors as 
mentioned above. 

The shape of the cross sections for both reactions shows 
wide structures peaked at about 1.5 GeV. Different in- 
termediate states contribute to the e+e~ — 47r cross 
sections. The observed bumps are sums of the contri- 
butions from the p(770), p(1450), and /3(1700) decays 
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FIG. 31 Comparison of the B ABAR res u lts o n the 
e^e^ — > 7r^7r~27r'' cross section (|Druzhininl. |2007|) with 
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into these intermediate states, which should be sepa- 
rated for a study of the excited p properties. Unfor- 
tunately, such analysis was performed at BABAR only 
at a qualitative level. The two- and three-pion mass dis- 
tributions for the 7r^7r~7r^7r~ final state are relatively 
well described by the model of the ai(1260)7r interme- 
diate state with a small contribution of the /o(1300)p 
state. This agrees with the ai tt dominanc e hypothesis 
suggested in Ref. ( Akhmetshin et al to describe 

the 47r dynamics at energies below 1.4 GeV. A strong de- 
viation from this hypothesis is observed in the 7r^7r~27r'^ 
channel. In addition to the expected ojtt^ and aiir con- 
tributions, a surprisingly large contribution of the 
intermediate state was observed. This is demonstrated 
in Fig. 1321 where the Att mass spectra for wtt, non-cj7r, 
and p~ intermediate states are shown together with 



the total mass spectrum for the e"* 



-27r>' 



action. The contributions of the different intermediate 
states were separated using simple conditions on Btt and 
27r invariant masses. It is seen that the p^ p^ cross sec- 
tion is more than a half of the non-wTr cross section at the 
energy about 1.7 GeV. For the 7r+7r~27r° masses higher 
than 2.5 GeV/c^ a clear signal of the /o(980) meson and 
a peak at the mass about 1.25 GeV/c^ (probably from 
the /2(1270) meson) are seen in the tt^tt" mass spectrum 
corresponding to the contributions of the /o(980)p and 
/2(1270)p intermediate states. 



F. e+e" ^ K+K-TT+ 

The e+e^ tt^ti^ reaction h as been never 

studie d before the BABA R experiment (|Aubert et all 
l2007bl: iLees et all 1201 1[ ). while the fully charged 
mod e was previously me asured with the DM1 detec- 
tor (jCordier et all Il982bl) but with an about 100 times 
smaller data set. The measured cross sections are shown 
in Figs. 1331 and [Ml The systematic uncertainties for these 
measurements are estimated to be at the (5-9)% level. 
The structures seen in the cross section energy depen- 
dence cannot be understood without analysis of interme- 
diate states involved. 

The distributions of the Kir invariant masses shown 
in Fig. I3S] indicate that the K*(9,mfK^iT^ and 
K* {9i92)^ K^-K^ (similar plots are not shown) intermedi- 
ate states dominate in these reactions. A small contribu- 
tion of the K^{IA3,Q)Kt: s taie is also s een (F ig.lgSTb')'). A 
special correlation study ( Lees et oil . 1201 1[ ) showed that 
the intermediate state with two K* , if * (892)^* (892), 
ii'*(892)i4'2*(1430), and K*{lAm)K*{lA?,Q), contributes 
less than 1% to the total reaction yield (the associ- 
ated X* (892)^1(1430) production is observed in J/V^ 
decays). Taking the numbers of events in the K* 
peaks for each cm. energy interval, the "inclusive" 
e+e- ^ K*{m2fK'K and e+e" -> K^{UmfKT: cross 
sections shown in Fig. [36] were extracted. Figures [3Wa) 
and ISTl b) show scatter plots of the reconstructed a) 
m(7r+7r~) and b) m(7r°7r°) versus m{K^K^) for se- 
lected events of the reactions e+e~ — )■ K'^K^tt^tt~ and 
e+e" — )■ K'^K^TT^TT^ , respectively. A clear 0(1020) sig- 
nal is seen in the K~^K~ invariant mass in both figures 
and is discussed in more detail below. A signal of the 
p{770) is observed in the tt+tt" invariant mass distribu- 
tion in Fig. I37r a). The tt+tt" invariant mass distribu- 
tion for K~Tr~^TT~ events not containing the K*{892) 
meson is shown in Fig. |3H] (a). The p(770) peak, prob- 
ably corresponding to the intermediate iiri(1230)^ and 
Ki{1400)^K^ states, is clearly seen in the tt+tt mass 
spectrum. In Fig. [35] (b) the "inclusive" cross section for 
the K^K^ p{770) reaction is presented. It is obtained 
by fitting the p(770) signal in the tt+tt^ invariant mass 
distributions for each cm. energy interval in Fig. 1381 

One o f the interesting ISR studies performed by 
BABA R (lAubert ef aLl.l2007bD and later reproduced by 
Belle ( Shen et all 12009 ) is extracting relatively small 
contributions of the (/)(1020)7r+7r-, 0(lO2O)7r°7r°, (</> -J> 
K^K^) intermediate states. Since the 0(1020) reso- 
nance is relatively narrow, the clean sample of (/itttt events 
can be easily separated. Figure [551 shows the m(7r+7r~) 
distribution for these events demonstrating a clear sig- 
nal from the /o(980) resonance and a bump at lower 
masses which can be interpreted as the /o(600) state. 
A similar plot is obtained for the tt^tt" invariant mass. 
These invariant mass distributions can be fitted with 
a superposition of two Breit-Wigner functions for the 
scalar /o(980) and /o(600) resonances as shown in Fig. [5^ 
The e+e" — )■ 0(lO2O)7r+7r~ cross section measured by 
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FIG. 32 The Att invariant mass spectrum for selected e+e" — >■ 7r+7r~27r° events (|Druzhininl l2007h (points with error bars) in 
comparison with the spectrum for non-aj7r''-events only (left) or with the spectrum for WTr^-events only (right). In the left plot 
the lowest histogram shows the contribution of the p~ intermediate state. 
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FIG. 33 The e^e~ — >■ K'^ K~-k^-k~ cross se ction measured 
with the BABAR detector (|Lees et a/.l . l201lh in comparison 
with t he only previous measurement by DM1 ijCordier et aZ.I . 
ll982bl V 



BABAR and Belle is shown in Fig. HOI Two resonance 
structures are seen at 1.7 GeV and at 2.1 GeV. The 
BABAR Collaboration investigated decay mechanisms 
for these structures and concluded that the second struc- 
ture decays only to the (/)(1020)/o(980) final state. The 
structure completely disappears if events associated with 
the /o(980) peak in the 77i(7r"'"7r~) distribution are re- 
moved. The first structure is associated with the (/)(1680), 
a radial excitation of the vector ss state. Its decays to 
<^(1020)/o(600) and 0(lO2O)/o(98O) are not forbidden. 

A simple VDM based model was suggested to de- 
scribe the observed e+e^ — > 0(lO2O)7r+7r~ cross section. 



FIG. 34 The e^e — ^ K 7r°7r° cross s ection measured 
with the BABAR detector (|Lees et aUbOllI '). 
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FIG. 35 (a) Scatter plots m{K~TT+) vs. m{K+n~), and (b) 
projection m{K^n^) plot (t wo entries per eve nt) for the re- 
action e'^e~ — >■ K~ n^n~ (|Lees et all 1201 il ). A sum over 
all accessible cm. energies is given. 
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FIG. 36 (a) The e+e" -» j^* ( 892)° j^7r, and (b) 
ifl (1430)''iir7r cross sections l|Lees ef a/.l . l20Hl ') obtained from 
the K*{S92f and 7^2(1430)° signals of Fig. IS^b), respec- 
tively. 
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FIG. 37 The scatter plots of the reconstructed a) m(7r^7r~) 
and b) m(7r''7r'') versus m{K^ K~) for selected events in the 
data. The vertical (horizont al) lines bound a (ji (/o(980)) sig- 
nal region l|Lees et all 1201 ll ) . . 

The model assumes that two vector mesons contribute 
to the cross section; one resonance associated with the 
(/)(1680) decays both to (?!)/o(600) and to (/)/o(980), while 
another referred to as F(2175) decays to 0/o(98O) only. 
Since the nominal 0(1680) mass lies below the 0/o(98O) 
threshold, the 0(1680) ^ 0/0(980) decay will reveal it- 
self as a smooth bump in the energy dependence of the 
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FIG. 39 The m{TT+ n~) distrib ution for the e+e" 
<?!>(1020)7r+7r" reaction (jLees et adhoil ). 
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FIG. 38 (a) The tt+tt" mass distribution for K+K-n+n- 
events {K* {892)Kn events are excluded); the solid curve rep- 
resents a fit using a signal Breit-Wigner function with p(770) 
parameters and a polynomial background (hatched area) . (b) 
The e^e~ — > K'^ K~ p{770) cross section obtained using the 
fitted nu mbers of p-meson events in each 25 MeV cm. energy 
interval ^^^di^Ol^. 



FIG. 40 The e'^e — >■ (pn'^n cross sections mea- 
sured with the BABA R (|Lees et all |2011| ) (circles) and 
Belle ijShen et all l2009l ) (squares) detectors. 



e+e~ -> 0/o(98O) cross section above 2 GeV. The result 
of the fit to the e+e^ 0(lO2O)7r"'"7r~ cross section with 
this model is shown in Fig. UT] It is clearly seen that 
the data above 2 GeV cannot be described with a con- 
tribution of the 0(1680) resonance only. An additional 
relatively narrow resonance y(2175) is needed to do this. 
The tongue below 2 GeV in the cross section for the re- 
action e+e~ 0(1680) 0/o(98O) in Fig.|4T]is due to 
the finite width of the /o(980) state. 

A relatively clean sample of 0/o(98O) events is se- 
lected using the requirement 0.85GeV/c^ < m(7r7r) < 1.1 
GeV/c^. The cross section for events of the K'tt^tt" 
mode, fitted with the model described above, is shown 
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FIG. 41 The fit to the e^e~ — >• (j)n^7v~ cross sec- 
tion ijLees et al\ . 1201 il l in the two-resonance model described 
in the text (soUd curve). The contribution of the first reso- 
nance ((^(1680)) is shown by the dashed line. The dotted line 
shows the first resonance contribution in the <^/()(980) decay 
mode only. 




FIG. 42 The cross section for e+e" -> (/)(1020)/o( 980) events 
selected with the cut 0.85 < m{mT) < 1.1 GeV/c^ (|Lees et all 
[2011). The solid curve is the result of the two-resonance fit. 
The dashed and dotted curves are the contributions of the 
(^(1680) 0/o(98O) and (?!)(1680) (^/o(600) decay channels, 
respectively. 



in Fig. 1321 The contribution of the y(2175) is seen 
much better with this selection. The comparison of the 
e+e~ — 7> (/)/o(980) cross sections measured by BABAR in 
two /o(980) decay modes, tt+tt" and tt^tt", is shown in 
Fig. SJl It is seen that two measurements agree. The 
fit of two modes gives the peak cross section, mass and 



FIG. 43 The e+e" (/>(1020)/o(980) cross section measured 
in the K~n^n~ (circles) a,nd K~ tt'^tt" (squares) final 
states by BABAR (|Lees et aUbOllf ). The solid and dashed 
curves represent the results of the two-resonance fit described 
in the text. 
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FIG. 44 The e+e~ - > ^(1020)/o(980) cross section measured 
by BeUe (jShen et al\ . l2009l ). The solid and dashed curves 
represent the results of the two-resonance fit described in the 
text. 



width of the resonance: 

ay = 0.104 ±0.025 nb, 
my = 2.179 ± 0.009 GeV/c^ 
Ty = 0.079 ±0.017 GeV. 



The 



— >■ (/) /n(980) cro s s sect ion measured in the 
Belle experiment (|Shen et all l2009l ) is shown in Fig. 
and also requires a resonance structure with similar pa- 
rameters. Some properties of this resonance, a relatively 
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FIG. 45 The e^e — >■ 2{K'^K ) cross section as a function 
of c m. energy measured with the BABAR detector using 
ISR l|Aubert et fflUl2007bl ). 



small width and absence of the (/)/o(600) decay, are un- 
usual. 

The nature of t his s tate is not 
clear ( Gomez- Avila et a/] . 120091 : iNapsuciale et all 
120071 ). One of the possible interpretations is that the 
y(2175) is a ssss four-quark state. Indeed, the /o(600) 
does not contain strange quarks, while the /o(980), 
strongly coupled with KK, definitely contains them. 



For the ssss state, the observed y(2175) 0/o(98O) is 
natural decay, while the not seen y(2175) — >■ 0/o(6OO) 
transition is suppressed by the OZI rule. The observation 
of the 1^(2175) decay to the 0?/ final state containing 
four s quarks, also supports this hypothesis. 



G. e+e" ^ 2{K+K~) 

The reaction e+e^ -> 2(K~^K~) was stu died for the 
first time by BABAR (|Aubert et a/.l . l2007bD . The mea- 
sured cross section is shown in Fig. 1451 The most sig- 
nificant structure in the cross sections is due to the J/V' 
decay. It is natural to expect that intermediate states 
for this reaction contain the </>(1020) meson which has a 
large decay rate to K^K~. Indeed, the strong meson 
peak is seen in the K^K^ invariant mass distribution 
shown in Fig. 1461 Since the (/) meson is present in al- 
most each four-kaon event, it is concluded that the re- 
action e+e^ 2{K^K^) is strongly dominated by the 
(/)K^K~ production. 

A study of events containing the meson was per- 
formed by BABAR. The K^K~ pair forming the me- 
son is selected by the requirement that its invariant mass 
is within ±10 MeV of the nominal mass. The in- 
variant mass distribution for the second pair is 
shown in Fig. HTT aV Figures ITTlb.c.d) show the cross 



FIG. 46 The K'^K invariant m ass distr i bution for selected 
e+e" 2{K'^K') events (Aube rt et oD. l2007bl ') (open his- 
togram, four entries per event), and that for the combina- 
tion in each event closest to the </)-meson mass (hatched his- 
togram) . 



section for events with the K~^K~ invariant mass in the 
regions 1, 2, and 3 indicated in Fig. ITfT a). An enhance- 
ment in the invariant mass spectrum near the 
K'^K~ threshold can be interpreted as being due to de- 
cay /o(980) — K^K~ . Therefore, the cross section for 
the region 1 is expected to have a structure similar to that 
observed in e+e~ — 0/o — K^K~mT (see Sec. IIII.Ft . 
The bump at 2.175 GeV is indeed seen in the cross sec- 
tion shown in Fig.lTflb). however, the data sample is too 
low to perform a quantitative analysis. 

The relatively narrow region 2 with 1.06 GeV/c^ < 
m{K^ K~) < 1.2 GeV/c^ is responsible for the spike 
seen at 2.25 GeV in Fig. |35l The spike is much more 
significant in Fig. |47] (c) showing the cross section for 
events from this mass region. There is no explanation of 
this structure. 

The peak in the K^K" mass spectrum near 1.5 
GeV/c^ is associated with the /2(1525). The region 3 
(1.45 GeV/c2 < m{K+K-) < 1.6 GeV/c^) ig chosen to 
select (f>f 2(1525) events. The cross section for this mass 
region is shown in Fig. ITTld) and exhibits a broad struc- 
ture at 2.7 GeV and a strong J/i/; signal. 



H. e+e — > 5 mesons 

The BABAR detector studied a number of 
ISR reactions with five hadrons in the final 

T^TT*', and 



state: 2(7r+7r ~)7r°, 2(7r+7r~)7?, K +K^tt+tt^ 
K+K-T:+n'r] (| Aubert al\ . \2007(h . 

The e'*"e~ — )■ 2(7r+7r~)7r° reaction has the largest 
cross section among the processes mentioned above. 
In the 7r"'"7r~7r'' invariant mass spectrum for this re- 
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FIG. 47 The K'^K invariant mass distribution for (f>K'^ K events (jAubert et al 1 l2007bl ) (a). Events from the J/ip -)> 
(j>K^K~ decay are excluded from the spectrum shown by the open histogram. The hatched histogram is for events from the 
J/tp decay. The numbered regions of the K^K~ mass spectrum are used to calculate the cross sections shown in the plots (b), 
(c), and (d) for the regions 1, 2, and 3, respectively. 
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action (Fig. |48)) clear signals of rj and uj mesons are 
seen corresponding to the wtt+tt" and r]TT~^Tr~ interme- 
diate states. The cross sections for t hese reactions w ere 
measured in direct e"'" e ~experimen t s (lAkhmetshin et 



Druzhinin et all 



2OOOI: lAntonelh aLl. Il988l Il992l: ICordier et all Il981 



1986( ). but BABAR data are signifi- 
cantly more accurate. The e+e" — ryTr+Tr" and e+e" — 
WTT+TT" cross sections measured by BABAR and in di- 
rect 6+ e~ experiments are shown in Fig. l49l and Fig. [50l 
respectively. 
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FIG. 48 The m(7 r +7r~7r° ) distribution for 2(7r+7r")7r° 
events l|Aubert et aZ.Ll2007d '). 
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FIG. 49 The e"'"e ~ — >■ T?7r^7 r ~ cros s section mea- 
sured by BABAR (|Aubert et all | 2007(J) in comparison 
with the direct e^e~measurements JA chmetshin et all I2OOOI : 
lAntonelh 6^0/! Il992l : iDruzhinin ef aU 119861 )^ 



+ 

s 

t 

+ 




Ee.„..(GeV) 



FIG. 50 The e^e — » uJiv'^T r cro ss section mea- 
sured by BABAR (jAubert aZ.I . l2007ct) m comparison 
with the direct e^e~measurcments ('Ak hmetshin et "al. I I2OOOI : 
lAntoneUi et oZ.I . [1988 : .Cordier et aL .198i[ k 




The two pions from the reaction e+e — ?■ 7771+ tt pre- 
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FIG. 51 The m(7r"'"7r ) distribution for selected uin'^n 
events in data (points with err or bars) and in simulation (his- 
togram) (|Aubert et am2007d '). 



dominantly form the p(770). In the two-pion invari- 
ant mass spectrum for the e^e~ — >■ u)tt~^7t~' reaction 
shown in Fig. [CTla clear /o(980) signal is observed. The 
contribution of the cl>/o(980) intermediate state was ex- 
tracted, and the cross section e+e~ — >■ a;/o(980) was mea- 
sured for the first time. It is shown in Fig. [52] The 
e+e^ — > WTT+TT^ cross section after subtraction of the 
a;/o(980) contribution is shown in Fig. [53l The cross 
section is fitted with a sum of of two resonances. The 
fit result is shown in Fig. [53] and listed in Table [iT] 
The obtained parameters are in good agreement with 
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FIG. 52 The e^e~ — » a;,fo(9 80) cross section measured by 
BABAR (jAubert et aZ.I . l2007d l. 
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FIG. 53 The fit with two Breit-Wigner functions to the 
Lon^ir~ cross section with t he aj/o(980) contribution sub- 
tracted (|Aubert et aZ.l . l2007d '). 

those obtained for the tt+tt^tt*' channel (see discussion 
in Sec HIO) . 

The WTT^TT" and rjir^TT^ intermediate states do not 
saturate the total e+e^ — > 2(7r+7r^)7r'' cross section as 
shown in Fig. 1541 The ?7i(7r+7r~) and m(7r^7r'^) distri- 
butions for e+e^ — >■ 2(7r"'"7r~)7r'' events with the wtt+tt^ 
and rjTT^TT^ contributions excluded are shown in Fig. 1551 
From the analysis of these two-pion mass distributions 
it was concluded that the dominant intermediate state 
for these events is p^p^n^. The pn mass spectrum also 
exhibits a resonance structure with the parameters: 

mx = 1.243 ± 0.012 ± 0.020 GeV/c^; 
Tx = 0.410 ±0.031 ±0.030 GeV. 




Km. (GeV) 

FIG. 5 4 The e+e" 2(7r+7r")7r° cross section (jAubert et all 
l2007d ) and contributions from cjtt^tt (squares) and 
rj-K'^-K~ (triangles). 
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FIG. 55 The m(7r^7r~) (points) and m(7r*7r°) (histogram) 
distributions for 2(7r^7r~)7r'' eve nts with the ujn^ 7v~ and 
rin^n~ contributions excluded (| Aubert et all [2007d l . 



The yield of the X(1240) state is consistent with 
the complete dominance of the quasi- two- body reaction 
e+e- ^ p(770)X(1240) ^ p^p'^Tr'f'. The best can- 
didates fo^^£l240) arejhe 7r(1300) or ai(1260) reso- 
nances ([Amsler et al\ . \200^ . 

The e+e" — > 2(7r+7r-)77 reaction was studied for the 
first time by BABAR. The measured cross section is 
shown in Fig. 1561 A rich internal structure is expected 
for the 47r77 final state. The four-pion mass distribution 
exhibits a wide resonance structure which can be a mix- 
ture of the known p(1450) and p(1700) resonances. Fig- 
ure [STIa) shows the r]Tr~^Tr~ mass distribution with two 
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FIG. 56 The e^e — » 2(-k^-k )rj cross section measured by 
BABAR (lAubert et flZ.I . l2007d ). 
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FIG. 57 (a) The in{r}'K^-K ) distribution for 2{-k^-k )r] events; 
(b) The e^e~ — ^ ?7(9 58)7r^7r~ cross sectio n and the result of 
the Breit-Wigner fit l|Aubert et a/.ll2007d '). 



narrow peaks. The lowest mass peak corresponds to the 

7r~ cross sec- 



77'(958). The measured e+t 



77'(958)7r-' 



tion is shown in Fig. I57r b). The resonance-like structure 
observed in the cross section energy dependence is fitted 
with a single Breit-Wigner function. The fitted resonance 
parameters are 

o-o 0.18 ± 0.07 nb, 
m^r = 1.99 ±0.08 GeV/c^ 
= 0.31 ±0.14 GeV. 

There is no entr y for these parani eters in the cur- 
rent PDG tables (jAmsler et all |2008| ). Taking into ac- 
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FIG. 58 (a) The m(r;7r+7r ) distribution for 2(7r"'"7r )r; events; 
(b) The e^e~ — >■ /i (12 85)71"*" tt" cross sectio n and the result 
of the Breit-Wigner fit (jAubert et aZ.I . l2007d ). 



Fig. [HHfb). The observed resonance structure has pa- 
rameters: 

(To = 1.00 ±0.18 ±0.15 nb, 

= 2.15 ± 0.04 ± 0.05 GeV/c^, 
= 0.35 ±0.04 ±0.05 GeV. 

The mass and width are close to those measured in the 
reaction e"'"e~ — > 77'(958)7r+7r^, but the cross section is 
significantly larger. The observed structure can also be 
assigned to the p(2150) resonance. 
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FIG. 59 Th e e e — » K K tt -k tt cross section measured 
by BABAR l|Aubert et oU 12007 j ). 



cross section for the reaction e+e ^• 
tt+tt^tt'^ shown in Fig. 1591 was also measured by 



count possible large systematic uncertainties on mass and 
width, the observed resonance can be interpret ed as the 
p(2150), extensively discussed in the past f Amsl er" et all 
I2008D . 

Another clear structure seen in the rjir^TT^ mass dis- 
tribution (Fig. ISTTa)) and shown in detail in Fig. ISSTa) 
was interpreted as the /i(1285) meson. The mea- 
sured e^e~ — > /i(1285)7r^7r~ cross section is shown in subprocesses. The former is in good agreement with that 



The 

K+R- 

BABAR for the first time. The three-pion and two- 
kaon invariant mass spectra for this reaction are shown in 
Fig. [SOI Clear rj and w signals are seen in the three-pion 
mass distribution and a strong (j) signal in the K~^K~ 
mass distribution. Figure [5T] shows the calculated cross 
sections for the e~^e~ — >■ (jyi] (a) and e~^e~ — )■ ujK^K~ (b) 
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FIG. 60 The m(7r+7r"7r") (a) a nd m(K^K-) (h) di stribu- 
tions for K+ tt+7v'tv° events l|Aubert et all l2007d '). The 
hatched histogram represents the estimated non-ISR back- 
ground. 
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FIG. 61 (a) The e^e — > (;6?7 cross section measur ed by 
BABAR in the /S'+A'"7r+7r"7r° (lAubert ef a/.L l2007d ') (cir- 
cles) and K'^K'yj (squares) (jAubert et all l2008bf l final 
states, (b) The cross sectio ns for the e^e~ — >■ uj K'^K~ pro- 
cess measured by BABAR ijAubert et aZ.Ll2007d '). 



obtained in the -> 77 mode ( Aubert aZ.I . [2008b[ ) . It 
is a first observation of the process e+e~ — > ujK^K~ . 
The reaction e+e" — >■ i'C"'" j^r~7r^ 7r~?7 w as also studied 



by BABAR in Ref. (jAubert et al.i . \2007(i ). The measured 
cross section is small and rises from threshold to a max- 
imum value of about 0.2 nb at 2.8 GeV, followed by a 
monotonic decrease with increasing energy. The clear sig- 
nal of the (prj' (958) intermediate state is observed in the 
K^K^ and tt+tt^t; mass distributions. Unfortunately, 
the (j )p'(958) invar i ant ma ss spectrum is not shown in 
Ref. (jAubert et aLl . l2007cl) . This spectrum is interesting 
since for the four-quark y(2175) resonance (Sec. IIII.Fp 
the decay to (l)r]'{958) is expected. 



I. 



3(7r+7r-), 2(7r+7r-7r°) 



and 



The reactions e+e^ — > 3(7r+7r 
2(7r+7r~7r°) were studied before in a number of 
experirn eni^s, but with limited data 



samples (iBacci aL . 198lt iBaldini et al . 1988 ; 



Bisello et ai . 19811 : Cosme et all . 1979t Esposito et al. 
198lh . The BABAR detector studied the six-pion 



production u sing the ISR method from the threshold 
to 4.5 GeV (jAubert et all I2006U) . As a resuh, the 



statistical and systematic uncertainties on the cross 
sections were dramatically reduced. 

An interesting feature of the 3(7r+7r^) final state is 
the presence, among many tt+tt^ combinations, of only 
one p(770)° per event. No other intermediate resonance 
signals were observed. For the 2(7r"'"7r~7r°) final state also 
one p(770) only, neutral or charged, per event is observed 
in the expected proportion 1:2. 

In the 2(7r+7r^7r°) final state, 7] and w signals are seen 
in the tt^tt^tt'^ invariant mass distribution. A small frac- 
tion of events corresponds to the associated production 
of the r] and w. Selecting these 77a; events the e~^e^ — > ujij 
cross section shown in Fig. |62| (left) was measured for the 
first time. The observed resonance structure which is ex- 
pected to be the a;(1650) is fitted with a Breit-Wigner 
function. The fitted curve is shown in Fig. |62] (left). The 
obtained resonance parameters are listed in Table jlTj to- 
gether with the resonance parameters obtained from the 
fits to e'^e^ — > Stt and e'^e" — >■ wtttt cross sections (see 
discussion in Sec. IIII.jp . Comparison of the wStt and 



2(7^^ 



uj{782)r] contributions with the total e+e 
cross section is shown in Fig. [62l right). 

The total e+e" 2(7r+7r"7r°) and e+e" ^ 3(7r+7r~) 
cross sections shown in Fig. |63| have very similar energy 
dependence. The ratio of the cross sections is almost 
constant over the energy range under study. Its average 
value is equal to 3.98 ± 0.06 ± 0.41. A dip structure just 
below 2 GeV in the six-p ion cross sect ion was observed 
in the DM2 experiment ( Baldini et all . [l988 ,) and then 
confirmed in the diffractive photoproduction of six pions 
in the FOCUS experiment (jFrabetti et a/.l . I2OOOI ) . Such a 
dip at 1.9 GeV was also observed in the total cross section 
of e"' "e~ annihilation into h adrons by the FENICE detec- 
tor (jAntonelli et all Il996i ). This structure in BABAR 
data is fitted using the Breit-Wigner function coherent 
with the smooth non-resonant background. The fitted 
curves for both cross sections are shown in Fig. 1631 The 
following "resonance" parameters are obtained: 



TOgTT = 1.88 ±0.03 GeV/c 

m4^27r" = 1.86 ± 0.02 GeV/c 

Fe^ = 0.13 ±0.03 GeV, 

^4^2-^0 = 0.16 ±0.02 GeV. 



The parameter values seem to be essentially indepen- 
dent of the final-state charge combination. These values 
may be also comp ared with tho se obtained in the FO- 
CUS experiment (Frabetti et a/.l . ,2000.) : m = 1.91 ±0.01 
GeV/c2, r = 0.037 ± 0.013 GeV. The mass values are 
consistent, but the widths obtained by BABAR are sub- 
stantially larger. Note that typical widths of known 
isovector resonances with mass near 2 GeV/c^ are 200- 
300 MeV. Since the obtained mass of the resonance 
structure is close to the double proton mass, it may be 
interp reted as a proton-ant i protq n subthreshold bound 
state (jPatta and O'Donnelll . [20031 ). 
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FIG. 62 The e — > wr; cross section (left) and c ontribution of e^e — >■ tJTr^Tr tt" (squares) and e^e — > ojiy (triangles) cross 
sections to all e+e" 2(7r+7r"7r°) events (right) ijAubert el oZJ . l2006bl ). 




Ec.™.(GeV) E^^(GeV) 

FIG. 63 The e^e~ — > 2(7r "'"7r~7r°) cross section (left) and the e+e~ — >■ 3(7r"''7r~) cross section (right) with a fit to the 
( ant i-) resonance function ijAubert et al l l2006bl ). A resonance responsible for the "dip" is shown shaded. 



J. Summary 

The BABAR ISR study covers the low energy range 
of e+e~ interactions from the 2tt threshold to 4.0-4.5 
GeV with exclusively measured cross sections for many 
processes. Figure [64] shows all exclusive cross sections 
measured by BABAR in a single plot. One can see that 
in most of the cases cross sections strongly depend on en- 
ergy and their central values vary by five orders of mag- 
nitude. 

One of the purposes of the BABAR ISR program was 



to measure the total hadronic cross section in the energy 
range below 2 GeV with improved accuracy ( DruzhiniDl . 



I2007D . To finalize this program, the cross sections at least 
for the 7r+7r-37rO, -k+tt-Att^, K+K', KsKl, KsKltttt, 
KsK^n^n^ final states should be additionally measured. 

Note that the total cross section value is not the di- 
rect sum of the cross sections shown in Fig. [Mj Each 
channel has internal subprocesses which include different 
resonances with different branching fractions to the ob- 
served final states. To perform a correct summation, each 
subchannel should be extracted separately and corrected 
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TABLE II Summary of the i:i;(1420)(tj') and ct;(1650)(a;") resonance parameters obtained from the fits described in the text. 
The values without errors were fixed in the fits. 



Fit 



Lur] (Aubert et al . 2006b l ujtt'^tt (Aubert et al. . 



3iT CAubert et al. 2004b'l PPG (Amsler et al. 20081 



aou,' (nb) 



w' f 



10^ 



r,e-B„'/(eV) 
m„/(GeV/c2) 
r„'(GeV) 
(j)^' (rad.) 
o-Q^" (nb) 



BeeB. 



w" f 



10° 



reeB„"/(eV) 

m„"(GeV/c2) 
r„"(GeV) 
4)^r, (rad.) 
rrow (nb) 
xVn.d.f. 



3.08±0.33 



1.645±0.008 
0.114±0.014 







1.01±0.29 
0.13±0.04 
17.5±5.4 
1.38±0.02±0.07 
0.13±0.05±0.01 

TV 

2.47±0.18 
0.47±0.04 
103.5±8.3 
1.667±0.013±0.006 
0.222±0.025±0.020 


102±67 

34.9/48 



0.82±0.08 

369 
1.350±0.030 
0.450±0.140 



1.3±0.2 
286 
1.660±0.011 
0.220±0.040 


PDG 



1.40-1.45 
0.180-0.250 



1.670±0.030 
0.315±0.035 



for the decay rate of an internal resonance. 

The exclusive ISR study of hadron production allows to 
investigate and improve our knowledge of excited states 
for light vector mesons. For most of them the parame- 
ters are still rather imprecise and new investigations are 
needed. 

For multihadron final states it may be difficult to iso- 
late the contributions of particular vector resonances due 
to presence of many interfering intermediate states. An 
example is the reaction e+e" tt 



tt^tt", to which the 



WTT , aiTT, p'^p~ intermediate states give dominant con- 
tributions. The two latter states contain wide resonances 
and strongly interfere. A partial- wave analysis is required 
to separate the sub-processes of the e+e" — >■ 7r"'"7r~7r°7r'' 
and e^e^ — > 7r+7r^7r+7r^ reactions. We hope that 
BABAR has enough data to perform such an analysis. 
This is necessary to separate contributions of two excited 
p states, p(1450) and p(1700), and determine their pa- 
rameters. A detailed study of intermediate mechanisms 
strongly benefits from the quasi-two-body character of 
some final states, but is much more difficult for multi- 
body final states. 

The BABAR ISR data on isoscalar channels already 
allow to improve parameters of excited (/> and co states. 
A global fit to the isovector and isoscalar components 
for the proces s e+e~ -» j^*(892 )J^, and the e+e" t/)?; 
cross section ( Buon et all . Il98^ was used to determine 
parameters of the (/)(1680) resonance (see Table[l|. 

In Table |TT] we summarize the results of the fits to 
the e^e^ — Stt, e~^e~ wtttt, and e~^e ~ — » ur] cross 
sections performed b y BABAR in Refs. ( Aubert et all . 
l2004bl l2006bl l20q7cD. and compare t hem with the cor- 
responding PDG (|Amsler et all I2008D parameters. A si- 
multaneous fit to all three channels could significantly 
improve the results and give additional information on 
relative decay rates. 

Due to numerous extensive studies of various exclusive 
cross sections, we have learned a lot about the total cross 
section of e~^e~ annihilation into hadrons and its com- 



ponents allowing a more precise estimation of hadronic 
vacuum polarization effects to be performed (see also Sec- 
tion |VlI|- 



IV. BARYON FORM FACTORS 
A. General formulae 



The cross section of the proces s e'^e — » BB, where B 
is a spin- 1/2 baryon, is given by ( Renardl . Il98ll) 



da_ 

an 



4s 



|Gm(s)P(1 + cos2 



1 



G£(s)psin2 



-|G£;(s)n sin^ 61 cos 2v? 

T 



(11) 

(12) 



where /3 = -y/l — Am^/s and mg are the baryon velocity 
{v/c) and mass, C = y/(l — e~^ ) with y = Tr amB/ 0\/s is 



the Coulomb correction factor ( TzaraL llQTOl) for charged 
baryons (G = 1 for neutral baryons) , t = m /4M| is the 
inverse helicity suppression factor, Gj\/ and Ge are the 
baryon magnetic and electric form factors. The number 
of form factors (two) corresponds to two BB states with 
different angular momenta: '^Si and ^Di. At the BB 
threshold the D-wa.ve state vanishes, and \Ge\ = \Gm\- 
At high y/s the terms containing Ge are suppressed by 
the helicity factor 1/t. With unpolarized beams the total 
cross section is 



a{s) 



Ana^PG 
3s 



\GM{s)\' + ^\GEi3r 
It 



(13) 



As discussed above (see Sections III.AI and III.Dl and 
Fig. [T]), the detection efficiency in the ISR measurement 
with a tagged photon has weak dependence on the an- 
gular distributions of final hadrons. In the case of the 
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FIG. 64 The cross sections of e^e — >■ hadrons measured with the BABAR detector via ISR. The results on e^e — >■ tt^tt tt' tt' 
are prehminary. 



dibaryon production this allows one to measure the to- 
tal cross section [Eq. ([T3| ] independently of the relation 
between the electric and magnetic form factors. The ra- 
tio of the form factors can then be determined from an 
analysis of the baryon angular distribution. In direct 
e+e" or pp experiments the range of the accessible po- 
lar angles is limited by the detector acceptance. In this 
case the cross section cannot be measured in a model- 
independent way. The detection efficiency is determined, 
and the proton magnetic form factor \Gm\ is extracted, 
usually under t he assumption that \G m\ ^\Ge\- In the 
BABAR paper (|Aubert et all l2006al ) on the ISR study 
of the reaction e^e~ — > pp the effective form factor is 
introduced as a linear combination of IGa/P and IGsP: 



\F{mf 



2t\Gm{s)\' + \Ge{s)? 
2t + 1 



(14) 



With the effective form factor the total cross section looks 
Hke 



(To (to) = 



1 



.(1+ )|F(m)|^ 



(15) 



Sto^ ^ ' 2t' 

The effective form factor defined in such a way allows an 
easy comparison of the results of the model-independent 
ISR measurement with \Gm\ obtained in direct e+e^ and 
pp experiments under the assumption that \Gm\ — \Ge\- 



The modulus of the ratio of the electric and mag- 
netic form factors can be determined from the analysis 
of the baryon polar angle distribution. This distribution 
can be presented as a sum of the terms proportional to 
|GAfP and IGgp. For the e+e" ppj cros s section the 
fully differential formula can be found in ( Czvz et all . 
I2004D . In this process the 9p dependences of the Ge 
and Gm terms are not strongly different from sin^ 9p and 
1 -I- cos^ 9p, describing the angular distributions for the 
electric and magnetic form factors in Eq. (|12l) . Note that 
in direct e^e^ experiments with transversely polarized 
beams a study of the proton azimuthal angle distribu- 
tion can improve Ge/Gm separation (see Eq. ([T2|)). 



A nonzero relative phase between the electric and 
magnetic form factors manifests itself in a polariza- 
tion of the outgoing baryons. In the reaction e+e^ — > 
BB this polarizatio n is p e rpe ndicular to the production 
plane (Dubnickov a et al\ . [l996) . For the ISR process 
e+e^ — >■ BB"i the polariza t ion observables are analyzed 
in Refs. (jCzvz e^ ad 120071: iKardapoltze^ |2007|) . In the 
case of the A A final state the A — >■ p7r~ decay can be 
used to measure the A polarization and hence the phase 
between the form factors. 
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B. Measurement of time-like baryon form factors 



2.1 GeV. 



Measurements of the e+e pp cro ss section have 
been p erformed in e~^ e~ experim ents (|Abhkim et al. 
20051: lAntonelh et ali [1998; .Bisello et al. 



IQQOl: ICastehano et oil . Il973l: iDelcourt et ali 



1983, 



1979; 



Pedlar ad l2005( ) with a (20-30)% precision. The 



cross section and the proton form factor were deduced 
assuming \Ge\ = \Gm\- More precise measurements 
of the proton form factor have bee n performed in 
pp — » e"'"e~ experiments jAmbro g iani et al\ . Il999l : 
[Armstrong et all Il993t iBardhTe^ alS.'^ m^. In the 
PS170 experiment (|Bardin et ali 11994 at LEAR, the 
proton form factor was measured from threshold {jpp 
annihilation at rest) up to a mass of 2.05 GeV/c^. 
The ratio \Ge/Gm\ was measured using the angular 
dependence of the cross section and was found to be 
compatible with unity. The LEAR data show a strong 
dependence of the form factor on pp mass near threshold, 
and a weak dependence in the range 1.95-2.05 GeV/c^ . 
The Fermilab e xperiments, E760 ([Armstrong et al\ . 
Il993l ) and E835 (jAmbrogiani et ali ll999|), show that 
the form factor decreases rapidly at higher masses, 
in agreement with the perturbative QCD prediction 
Gm oc alirn^) / m,'^ . 

Experimental information on the reactions e^e~ — 
AA, TPyP, AS" is very scarce. The e+e~ AA cross 
section is measured to be 100^35 pb at 2.386 GeV, and 
at the same energy the upper limits for e'^e'^ — ?> 'TPYP 
(< 120 pb) and e+e ~ -» AS° i<J^ P b) cross sections 
have been obtained ( Bisello et aLl . ll99"of) . 



C. e+e — s> 



The first ISR baryon experiment was the measure- 
ment of the proton-antip roton production cross sec- 
tion (|Aubert et ali [2006a|) by BABAR. The measured 
e+e" — )■ pp cross section shown in Figs. IHS] and [51] is 
almost flat near the threshold, and then decreases from 1 
nb to about 1 pb at 4.5 GeV. There are two rapid drops 
of the cross section near 2.15 and 2.9 GeV. The BABAR 
proton form factor data presented in Fig. [HTl in gen- 
eral, agree with the previous measurements. Figure [68] 
shows an expanded view of the near-threshold region. 
The BABAR measuremen t confirms the PS 170 observa- 
tion ( Bardin et ali [l994[) of the significant increase in 
the form factor for energies approaching the pp thresh- 
old. The proton form factor reaches about 0.6 at the 
threshold. 

A study of the proton angular distribution allows one 
to extract the value of the ratio of the electric and mag- 
netic form factors \Ge/Gm\- The results of the BABAR 
\Ge/Gm\ measurement are shown in Fi g. [Ml in com- 
parison with the data obtained at LEAR ( Bardin et ali 
]T994). In disagreement with the LEAR result, the 
BABAR data indicate that \Ge/Gm\ significantly ex- 
ceeds unity in the energy range between threshold and 
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FIG. 65 The e e — >■ pp cross section measured by 
BABAR ([Aubert et d.L [2006ai ) . 
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FIG. 66 The e'^e — » pp cross section measured by 
BAB AR (lAubert et ali l2006al ') in comparison with the BES 
data ([Ablikim et ali [2005] ) 



D. e+e" 



AA7 



-> AA cross sect i on mea sured by the 
BABAR detector (jAubert et ali [2007d[ ) is shown in 
Fig |70l in comparison wit h the only previous measure- 
ment (|Bisello et ali [1990^ . The BABAR measurement 
is based on about 200 AA events selected in the decay 
mode A — >■ p-K. 

The measured A effective form factor is shown in 
FiglTT] The ratio \Ge/Gm\ is found to be consistent 
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FIG. 67 Th e proton form fact or measured in d if ferent 

experiments ([Ablikirn et all l2005l: lAm brogiani et al[ 1 19991 : 



1993; Aube rt et al. 
1983; Dclcourt et al. 



AntoncUi et al, 1998; Armstrong et 
2006a; Bardin et al, 1994; BiscUo et al, 
19791 : iPedlar ef all 12005'). The sohd hne represents the QCD 
fit described in the text. 
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FIG. 68 The proton 

region (lAntonelli et ali . . 

iBardin ef aZ.I . Il994l : iDelcourt et aZLIl979l ) 



form fact or in the near -th reshold 
19981: lAubert et all l2006al : 



with unity. The use of the A — >■ pvr decay allows to mea- 
sure the relative phase 0a between the complex Ge and 
Gm form factors. A non-zero 4>\ leads to polarization 
C of the outgoing baryons. The value of C is extracted 
from the analysis of the proton angular distribution in 
the A — 7> pTT decay. The measured cos^p^ distribution, 
where 9p(^ is the angle between the A polarization vector 
and the proton momentum in the A rest frame, is shown 
in Fig l72l No cos 9pQ distribution asymmetry correspond- 
ing to the non-zero polarization is seen. Because of the 



FIG. 69 The proton \G e/Gm\ ratio measured by 
BABAR (lAubert et aZ.!. 1200681) f blacl<: points) compared with 
LEAR data (|Bardin et all\l994 ) (open circles). 



limited data sample only a very weak limit on the phase 
between the Ge and Gm has been set for the A hyperon: 
-0.76 < sin0A < 0.98. 
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FIG. 70 The e'^e — >■ AA cross section measured by 
BABAR (lAubert et all l2007dt) in comparison with the DM2 
measurement (jBisello et all Il990h . 



e e — > 



E°E°, AE° (E°A) 



The BABAR measuremei it of the S° and S°A form 
factors is described in Ref. ( Aubert et all [2007df ) . The 
decay chain E° — > A7 — > pirj is used to reconstruct 
About 20 candidate events were selected for each 
of the ISR reactions, e+e~ I]°i;"7 and e+e~ -> I^^Aj. 
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FIG. 71 The baryon fo rm fectors measured by 

BABAR (|Aubert et all l2006al . 12007(31 ) versus the dibaryon 



invariant mass. 
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FIG. 72 The cosd„c distribution in the e+e 



AA pro- 





cess (jAubert et aZ.I . l2007dr ). The line represents the result of 
the fit to data with a first-order polynomial. 



The effective TP and S°A form factors are shown in 
FigEl The corresponding values of the e+e" — )■ S°E" 
and e+e~ — t- E^A cross sections are about 40 pb near the 
reaction thresholds. It is seen that the A, E° and E"A 
form factors are of the same order. 



F. Summary 

The baryon form factors are a su bject of vari- 
ous p henomenological models (see Ref. ( Baldini et all 
|2009[ ) and references therein). QCD predicts for the 
baryon form factor the asymptotic behavior F{q^) ~ 



"^((7^ )/<?■* (|Chernvak and Zhitnitskvl . Il977l ). Compari- 
son of this prediction with the data on the proton form 
factor is shown in Fig. [67l It is seen that the asymptotic 
regime is reached at energies higher than 3 GeV. 

The remarkable feature of the process e~''e^ — ?> pp 
is a nearly flat cross section in the 200-MeV region 
abov e the pp thre s hold. This feature is explained in 
Ref. ( Baldini et al\ . |2009| ) by the opposite trends in the 
energy dependence of the 5- wave and D-wave contribu- 
tions. 

A natural explanation for the sharp increase of the 
proton form factor in the vicinity of the pp threshold 
is the final state interactio n of the proton and antipro - 
ton (see, for example, Ref. ( Dmitriev and Milsteinl . I2OO7I ) 
and references therein). Another possibility is a con- 
tribution of the vector-meson state located just below 
the pp thresh old. This stat e is observed in the reaction 
Gtt (lAubert et all l2006bt iBaldini et all . Il988l: 
|2000D . 



iFrabetti et al. 



The rapid drop of the cross section at 2.15 GeV may 
be a manifestation of the isovector state p(2150), which 
is seen in the r eactions e"'"e~ — > 7?V "'"7r~ and e"'"e~ — > 
/i(1285)7r+7r- (jAubert et all l2007cl ). The drop in the 
cross section near 2.9 GeV is still not understood. 



The 



— 7> AA cross section has some features simi- 



lar to those for the process e+e^ — > pp. In the energy re- 
gion of about 200 MeV above threshold the AA cross sec- 
tion is flat; the Ge/Gm ratio is consistent with that mea- 
sured by BABAR for the e"''e~ — )• pp. The attempt to ex- 
plain the unusual energy depende nce of the e"'"e~ — » A A 



cross section was made in Ref. ( Baldini et al ], l2009t ). 

A fit to the A form factor with the power-law func- 
tion const/q"- (FigESl) gives n = 9.2 ± 0.3. This n 
value strongly differs from the QCD asymptotic predic- 
tion 71 — 4. Similarly to the pp case, the asymptotic 
regime is not reached for e+e" — >■ AA at the energies 
below 3 GeV. 

The cross sections of the processes e+e" — )■ YPyP 
and e+e~ — )■ AE'^(S°A) have been measured with large 
errors. The corresponding and E°A form factors 
fFig |7T|) show a monotonic decrease starting just from 
the threshold. A fit to form factor data with the power- 
law function gives n > 4 but with large errors. 

It is interesting to compare the measured form factors 
with each other and with the QCD prediction for th e 
asymptotic form factor ratios ( Chernvak et all [1989'): 
Fp = 4.1fA, i^so = -1.18Fa, Feoa = -2.34Ja. From 
comparison of the form factors in FiglTTl it is seen that 
the prediction works (possibly accidentally) only for the 
ratio of the A and E'^ form factors. The ratio Fa/Fp falls 
with energy. In the highest energy interval 2.8-3.00 GeV 
the ratio is equal to 0.3+q'3 and agrees with the asymp- 
totic value 0.24. This is an indication that the asymptotic 
regime is reached just above 3 GeV. 

The BABAR experiment shows that the ISR method 
is well suited for the measurement of baryon form fac- 
tors. Future Super _B-factories as well as already running 
BEPC e+e~ collider will make possible measurements 
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particle identification in particular, are much better com- 
pared to the BESII detector. As a result, the detector 
efficiency and the integrated luminosity are determined 
with lower systematic errors. A typical systematic uncer- 
tainty of the BABAR measurement is 3-5%, while BESII 
usually quotes 10-15%. The lower systematic error makes 
ISR results on many J/-0 decays competitive with BESII 
and other previous measurements. Practically all decays 
with the rates about 10"'^ and higher can be measured 
via ISR with better overall accuracy. Moreover, because 
of excellent particle identification, many J/ij} and ip{2S) 
decays with kaons in the final state have been studied 
using ISR for the first time. 

In the BABAR experiment the ISR method enabled 
to measure a few tens of J/ijj and ip{2S) decays with the 
best-to-date accuracy and discover about 20 new decays 
of these resonances. 



FIG. 73 A fit of the A form factor (jAubert et aU 12007^ 1 with 
the power-law function F ~ const /q" and with the QCD- 
inspired function F ~ const /q'^. 



of the form factors, especially for the proton, with un- 
precedented accuracy. High-precision measurements of 
the proton form factor are also planned in the PANDA 
experiment. . 



V. DECAYS OF THE J/^p AND ^^{23) 

For all the processes described in the previous sections, 
clear narrow peaks are seen in the energy dependence of 
the cross sections corresponding to the J/tp and il^{2S) 
decays. The Born cross section for the ISR production 
of a narrow resonance, for example, the J / ip, decaying to 
the fin al state h can be calculated using ()Aubert et al\ . 
l2004aD 



-)B{J/i>^h) 



(16) 

where mj/^ and Y{J jil) — >• e+e ) are the mass and elec- 
tronic width of the J/V' meson, xj/^ = 1 



/s, and 

B{J /il) — > h) is the branching fraction for the J /ip de- 
cay to the final state h. The function Wo is described in 
Sec. lI.CI bv Eq.©. Therefore, a measurement of the num- 
ber oi J/ip ^ h decays in the ISR process e+e^ 
determines the product of the electronic width and the 
branching fraction: r{J/tp e'^e~)B{J /ip — h). 

The total cross section for the process e+e~ — 7J/V' 
with a tagged ISR photon (^o = 30°) is about 3.4 pb. 
With the integrated luminosity of ^ 500 fb^^ collected 
by the BABAR detector it corresponds to about 1.7 mil- 
lion produced J/i/j's. This number is significantly smaller 
than, for example, about 60 million J/tp^s produced in 
the BESII experiment at the BEPC e+e~ collider. How- 
ever, the general quality of the BABAR detector and its 



A. Leptonic decays 



The BABAR (lAubert all l2004a^ and 
CLEO (lAdams et al\ . I2006D collaborations per- 
formed a study of the J ftp production in the reaction 
e+e^ — >■ 11'^ fj,^^. The dimuon mass spectrum for this 
reaction obtained by BABAR is shown in Fig. [TH The 
signal of the ISR J/ip production is well seen in the 
mass spectrum. The nonresonant spectrum is due to 
muon pair production in the process e~^e~ — ^^fj,~j, 
where photon can be emitted by both initial electrons 
and final muons. Since the dimuon decay of the J/tp 
meson proceeds through a single photon transition, 
J/ip — >■ 7* — > fX^ ^~ , the angular and momentum dis- 
tributions for events from the J/ip peak are completely 
identical to those for the ISR part of nonresonant events. 
The idea of the BABAR measurement is to determine 
the ratio of the number of J/ ip events to the level of the 
nonresonant spectrum which is well known theoretically. 
The dimuon spectrum of Fig. [73] has been fit with a 
function taking into account the energy dependence of 
the nonresonant cross section and the experimental J/ip 
line shape. The ratio 



TV 



^ • Am 

dm 



(17) 



was the main fit parameter. After substituting cross sec- 
tions for the numbers of events, this ratio can be rewrit- 
ten 



J/ ■4' 



do-" 



dm 



• Af 



d<tai/dm 
daitn/dm 



(18) 



Detector acceptances and ISR radiative corrections, 
which are the same for the nonresonant ISR and J/ip 
contributions to the reaction e+e^ — > /i+/x^7, cancel in 
the ratio. The total nonresonant cross section includes 
the FSR contribution, which is parameterized in terms 
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of K, the ratio of the visible nonresonant total and ISR- 
only (FSR switched off) cross sections. Since BABAR 
selects events with the photon emitted at a large angle, 
the FSR contribution is relatively large. Using simulated 
events, the coefficient K = 1.11 ± 0.01 (statistical error 
only) is determined for the selection criteria used. 



Si 000 



by CLEO at the il'{2S) resonance were used to select 



500 




3.2 



3.4 

M^^ (GeV/c') 



FIG. 74 The Ai^/x mass spectrum in the J/il) reKion for 
selecte d events of the reaction e~^e~ — > pT^ (jAubert et all 
l2004al ). The curve is the result of the fit described in the text. 



a clean sample of '4>{2S) — >■ J/i/jTr+Tr , J/i/j 



events. The CLEO result for the product of the J/'0 
parameters is 

Tee • = 0.3384 ± 0.0058 ± 0.0071 keV. 

Despite the analysis improvements described above, the 
systematic error of the CLEO result was not reduced 
compared to the BABAR measurement. However, the 
sources of the systematic uncertainties are different for 
the two measurements. So the results can be considered 
as completely independent. 

The data collected at xfs = 3.77 GeV were used 



by C LEO to study ■0(25) ISR production ()Adam et all 
I2006D . ISR il){2S) events were selected in the decay 
modes ipi^S) — > n^TT^J/tp, -K^ir^J/tp, and -qJ/ip with 
the J/ijj decaying to the lepton pair, e+e" or //"""/x". 
From the number of the V'(2S') events the products 
T{il}{2S) e+e-)B{ip{2S) XJ/^), where X = tt+tt", 
tt'^tt'^, and rj, were obtained. Since the branching frac- 
tions for these decay modes are known with 1.5-2% ac- 
curacy (|Amsler et all 120081) . the measurement of the 
products can by used to improve accuracy of the ^p(2S) 
electronic width. The CLEO result dominates in the 
curre nt PPG value TM 2S) e+e") = 2.36 ± 0.04 
keV (jAmsler et a^.l . |2008[) . 



The result of the fit is shown in Fig. [74l The value 
r = 18.94 ± 0.44 is found with x^/ndf = 122/144. From 
the product r ■ K ^ 21.03 ± 0.49 ± 0.47 the cross section 
aji^ = 2124 ± 49 ± 47 fb and the product of the J/^ 
parameters 



Bf,f, = 0.3301 ± 0.0077 ± 0.0073 keV 



are determined. The main sources of the systematic error 
quoted are uncertainties in the J/'0 line shape and the 
coefficient K , both due to imperfect simulation of the 
detector response. 

Us ing the values for Bee and B^^ ( Eidelman et al\ . 
12004), which are w ell measured in t he cascade ip{2S) — ?> 
J/t/^TT+TT^ decays ()Bai et all Il998() . the electronic and 
total widths of the J/ip meson were derived, 

Fee = 5.61 ± 0.20 keV, F = 94.7 ± 4.4 keV . 

These were the best-to-date measurements of the J/ip 
parameters. 

The BABAR measurement was improved by CLEO. 
With the integrated luminocity of 281 pb~^ collected at 
y/s = 3.77 GeV about 13 x 10^ ISR produced J/ip -> 
^J^' li~ events were selected (compared to 8 x 10'^ J/?/; — 
events in the BABAR measurement based on a 
88 fb~^ data sample). Since CLEO used the untagged 
approach, the FSR contribution to the nonresonant cross 
section was significantly reduced. The second important 
improvement of the method was that the J line shape 
was extracted from data. To do this, the data collected 



B. Decays to light mesons and baryons 

A systematic study of the J/i/' and i{j{2S) decays 
to light hadrons was performed in t he BA B AR exp er- 
iment (|Aubert et all l2004bl . l2005bl . l2006allbl . l2007cl ldh. 



An example of the J/ijj signal for J/ip ^ tt^tt tt , one 
of the most probable J /ij} decay modes, is shown in 
Fig. [75] (lAubert et all l2004b^ . It is seen that the nonres- 
onant background is small. From the number of events 
at the peak, the product V{J/ip e'^e~')B{J/tp 
Stt) = 0.122 ± 0.005 ± 0.008 keV was determined. Using 
the J/ip electronic width value, known from the ISR 
study of the J/4' — ^ decay, the branching fraction 

B{J/iIj Stt) = (2.18±0.19)% was calculated, which dif- 
fered by about 50% from the PPG value, (1 . 47 ±0. 13)%, 
available when the analysis (jAubert et al. I, l2004bD was 
carried out. Similar deviation was observed in the BES 
experiment (|Bai et all \2004 l where B{J/tp Sn) = 
(2.10 ± 0.11)% was obtained. 

Another example of a J/i/; decay mode with a rather 
high probability, which was studied using ISR, is shown 



[n'TT TT 



in Fig. [TBI where the signals of J/ip 
and '0 (25') — > 2(7r+7r^)7r'' are clearly seen (jAubert et all 
l2Q07d) . A gain, by determining the number of peak events 
over the nonresonant background and using Eq. (|16p . the 
product T{J/4' — > e+e~)B{J/ip — > Tr+vr^Tr+vr^Tr'') = 
(3.03 ± 0.05 ± 0.18) X 10"'' keV was determined. The 
value of the J/iIj — !■ 7r+7r^7r+7r~7r'^ branching fraction ob- 
tained from this product, (5.46±0.09±0.34)%, differed by 
about 5a from the PPG value, (3.37 ± 0.26)%, available 
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TABL E III Measurements of the J/ip and ^{2S) branc hing fractions via I SRat BABAR ijAubert et aZ.l . l2004bl . l2005bl . l2006al lbl. 
l2007d ldl) compared to the current world-average values (|Nakamura et aZl . |2010| ). 





Measured 


Measured 


J/i/) or tp{2S) branching 


fraction (10"^) 




quantity 


value (eV) 


BABAR 


PDG-2010 



^ ee ^J/^—)--K+7\ 
r ee O j 



- ee 

- ee 
ee 

- ee 

- ee 



-'j/i/>->2(ir+7r-) 
J/i/i->2(7r+7r-)7r0 
J/i/'->3(7r+7r-) 

J/V->-Ff+X-7r+ir- 
J/i/>->if+if "7r°7r0 
J/ V -!■ if + -ff " -ff + if ^ 
J/ip^K+K~Tr+n~n° 
J/-4>^K+K-2(w+n-) 



-pJ/tpfO f2 

Tee ^J/tp^-qir + ir- ^rj^Z-n 



ee 



B 



Tee Sj/i/)->0?7^,^->K+if- ^iJ^Stt 
-Lee Ojj^^jf^n- 



.27 



J- ee 'Jj/^ 
Tee Sj/^^if+K- 7r+ir-r;S»;->77 



^ee B j^^^i^jf,0x*0+C.C.)^K*0^K7rBji*0^KTr 



1 ee 



- ee 



J/!/'->(ifOif*(892)0+C.C.) 
J/!/'^(if+if*(892)-+C.C.) 



i ee '~> 



^<I>^K+K- 



Tee Sj/',j,^4,K+ K-^<t,^K+ K- 



B 



ee 

■pJ/i'r? 
^ ee ^ 



Tee Bj/^^^fgB^^ji+x-Bfg^^O^O 



)^il>^K+ K- 



122 ± 5 ± 8 
19.5 ± 1.4 ± 1.3 
303 ± 5 ± 18 

23.7 ± 1.6 ± 1.4 
89 ± 5 ± 10 

37.9 ± 0.8 ± 1.1 

11.8 ±0.8 ±0.9 
4.00 ±0.33 ±0.29 

107±4±6 

27.5 ± 2.3 ± 1.7 
47.8 ±3.1 ±3.2 

22±3±2 
0.51 ±0.22 ±0.03 
5.16 ±0.85 ±0.39 
0.84 ±0.37 ±0.05 
4.8 ±0.7 ±0.3 
3.3 ± 1.3 ±0.2 
10.2 ±1.3 ±0.8 
8.59 ±0.36 ±0.27 

26.6 ±2.5 ±1.5 
29.0 ± 1.7 ± 1.3 

0.57 ±0.15 ±0.03 
2.19 ±0.23 ±0.07 
1.36 ±0.27 ±0.07 
2.26 ±0.26 ±0.16 
0.69 ±0.11 ±0.05 
0.48 ±0.12 ±0.05 
4.7 ±0.9 ±0.9 



21.8 ± 1.0 ± 1.6 
3.70 ±0.26 ±0.37 

54.6 ±0.9 ±3.4 
4.40 ±0.29 ±0.29 

16.5 ± 1.0 ± 1.8 
6.84 ±0.15 ±0.27 
2.12 ±0.15 ±0.18 
0.72 ±0.06 ±0.05 

19.2 ± 0.8 ± 1.5 
5.09 ±0.42 ±0.35 

9.7 ±0.6 ±0.6 

4.1 ±0.6 ±0.4 
0.40 ±0.17 ±0.03 

2.35 ±0.39 ±0.20 
1.4 ±0.6 ±0.1 

0.87 ±0.13 ±0.07 

1.36 ±0.50 ±0.10 

4.7 ±0.6 ±0.4 
6.98 ±0.29 ±0.21 

4.8 ±0.5 ±0.3 

5.2 ±0.3 ±0.2 
0.23 ±0.06 ±0.01 
0.81 ±0.08 ±0.03 
0.50 ±0.10 ±0.03 
1.67 ±0.19 ±0.12' 
0.38 ± 0.06 ± 0.02 
0.53 ±0.13 ±0.05 
1.77 ±0.35 ±0.12 



20.7± 1.2(S = 1.2)" 
3.55 ±0.23 
41 ±5(5 = 2.4) 
4.3 ±0.4 
16.2 ±2.1 
6.6 ±0.5 
2.45 ±0.32 
0.76 ±0.09 
17.9 ±2.9(5" = 2.2) 
4.7±0.7(S = 1.3) 
8.6 ±0.7(5 = 1.1) 
4.0 ±0.7 
0.40 ±0.17 
2.29 ±0.24 
0.75 ±0.08(5 = 1.5) 
0.87 ±0.15 
1.36 ±0.51 
4.67 ± 0.70 

6.0 ±0.6 
4.39 ±0.31 
5.12 ±0.30 
0.23 ±0.07 
0.94 ±0.09(5 = 1.2) 
0.56 ±0.16 
1.83 ±0.24 
0.32 ±0.09(5 = 1.9) 
0.32 ±0.09(5 = 1.9) 
1.66 ±0.23 



2S)-+2(7r+ir-)7rO 

2S)-»2(7r+ir-7rC) 

2S)-»if +if-2(ir+7r-) 

2S)^J/ilnr+Tr- ^J/ip^ST: 

2S)^-cJ7r+7r- S^^-Svr 

2S)-»J/Vl7^'7^37r^3j/^_»^+ 

2S)->2(7r+7r-)?7^»7->77 

2S)^K+K-TV+TV-TV° 

2S)^-if+ if- 71+71-17^17^77 

2S)^-if+if-7r+7r- 

2S)^-if+if-7r07rO 

2S)^-if+if-if+if- 



2S)-+(^>7r+7i 

.B, 



B, 



Bt 



2S)-^4>fo'~'4>^K+ K- 1-' fa^TT 



29.7 ±2.2 ±1.8 
11.2 ±3.3 ±1.3 
4.4 ±2.1 ±0.3 
18.6 ±1.2 ±1.1 
2.69 ±0.73 ±0.16 
1.11 ±0.33 ±0.07 
1.13 ±0.55 ±0.08 
4.4 ± 1.3 ±0.3 
1.2 ±0.7 ±0.1 
1.92 ±0.30 ±0.06 
0.60 ±0.31 ±0.03 
0.22 ±0.10 ±0.02 
0.27 ± 0.09 ± 0.02 
0.17 ±0.06 ±0.02 



12.0 ±0.9 ±0.7 
5.3 ± 1.6 ±0.6 

2.1 ± 1.0 ±0.2 
23.6 ± 1.6 ± 1.6 

1.22 ±0.33 ±0.07 
33.4 ±9.9 ±2.0 

1.2 ±0.6 ±0.1 
1.8 ±0.5 ±0.1 

1.3 ±0.7 ±0.1 
0.81 ±0.13 ±0.03 
0.25 ±0.13 ±0.02 
0.09 ± 0.04 ± 0.01 
0.35 ±0.12 ±0.01 
0.22 ±0.08 ±0.02 



2.9 ± 1.0(5 = 4.6) 
5.3 ± 1.7 
1.9 ±0.9 
20.7 ± 1.2(5 = 1.2)" 
0.73 ±0.12 
32.8 ±0.7 

1.2 ±0.6 
1.26 ±0.09(5 = 1.2) 

1.3 ±0.7 
0.75 ±0.09(5 = 1.9) 

0.25 ±0.13 
0.060 ±0.014 
0.117 ±0.029(5 = 1.7) 
0.068 ± 0.024(5 = 1.1) 



is a PDG scale factor. 



J/tlj^4>KK 
^J/V'->37r- 
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or 777r+7r final states. It was found to be less than 3%. 
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FIG. 77 (a) The three-pion combination closest to J/^ mass 
versus five-pion mass, (b) The five-pion mass for events with 
the t hree-pion mas s in th e ±50 MeV window around J/tp 
mass (|Aubert et fflm2007d ). 



FIG. 75 The Stt mass spectrum for selected e^e~ — > 
7r'^7r~7 r°7 data events in the vicinity of the J/^p reso- 
nance (|Aubert et al. 
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FIG. 
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mass distribution for ISR 



76 The 

e^e~ — » 2(7r"'"7r~)7r°7 ev ents in the J/ip-4>{2S) mass re- 
gion (jAubert et all l2007d ') . 



when the analysis ( Aubert et g/J . [20073 ) was carried out. 
As was shown in Sec. IIII.Hi the five-pion final state in- 
cludes production of many intermediate resonances. All 
of them are seen in the J/ijj — Stt decay. This may be 
a source of the systematic error unaccounted in previous 
measurements of the decay. The detection efficiency in 
the ISR method with a tagged photon is weakly sensi- 
tive to the dynamics of the J/tp — >■ Stt decay. The model 
uncertainty in the detection efficie ncy for the BABAR 
measurement ( Aubert et al\ . \2007(i) was estimeted from 
the difference in efficiency values for phase-space gener- 
ated five-pion events and events generated for the wtt+tt" 



A part of events from the ip{2S) peak comes from the 
decay chain ip{2S) — >■ J/ipTr'^7T~ 2(7r"'"7r~)7r° with the 
J/ijj decaying to three pions. To select these events, the 
tt+tt^tt'^ combination with the invariant mass closest to 
the J/ip mass is chosen. Figure [TZT a) shows the scatter 
plot of this three-pion mass versus the five-pion mass. A 
clear signal from the above decay chain is seen. The five- 
pion mass spectrum for events with the 7r"'"7r~7r'^ mass 
within the ±0.05 GeV/c^ window around the J/^p mass 
is shown in Fig. iTTl b). From the fit to the mass spec- 
trum with a double-Gaussian function the number of de- 
tected ^p{2S) J/^n'^n~ — 7^ 2(7r+7r^)7r" events was de- 
termined to be 256 ± 17, and the triple product 



Bitp{2S) J/tpTT+TT-)B{J/tlj tt+tt-tt") 
(1.86 ± 0.12 ± 0.11) X IQ-^ keV 



e+e") 



was obtained. By using the world-avera ge r(^(25') — > ■ 
e^_e~) and B{ip{2S) — J/V'tt+tt") values ( Amsler et all . 
l2008h . the branching fraction S(J/-0 ^ tt+tt-tt") = 
(2.36±0.16±0.16)% was obtained, which is in good agree- 
ment with the BABAR measurement in the Stt final state : 
BjJ/'tp tt+tt-ttO) = (2.18 ± 0.19)% (lAubert et al\ . 
l2004bD . This, in particular, confirms the correctness of 
the normalization procedure used for the measurement 
of S(J/V' 57r). 

Table Hill presents measurements of the J/tp and tp{2S) 
decay rates performed with the BABAR detector via ISR 
for many multihadron final states. The current PPG 
value s are shown in the last column ( Nakamura et all . 
I2010D . In most of the cases these values are close to 
those of BABAR emphasizing their importance. Note 
also that in a few cases the scale factor is significantly 
higher than one indicating a large difference between the 
BABAR measurement and previous results. 

As can be seen from Table IIIIl the J/ tp decay rates 
to even numbers of pions (47r, Gtt . . . ) are much smaller 
compared to the decays to odd numbers of pions. Indeed, 
a strong decay of the J/tp to an even number of pions is 
forbidden by G-parity conservation. It is expected that 
this decay is dominated by a single photon transition. 
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FIG. 78 The 2(7r+7r-) mass distribution for ISR-produced 
e^e^ — > 2(7 r^7r~) events i n the m ass region around the J /ip 
and ip{2S) (jAubert et all l2005bl ): there are clear signals at 
the J ftp and ip{2S) mass positions. The latter is dominatesby 
the ip{2S) — > J/il)TT^TT~ — > fi'^ fj,~ -K^TT~ transition; selected 
events with two muons from the J/ip decays are shown by the 
shaded histogram. 



J/^p — > 7* — > riTT. No such suppression occurs for the 
strong J/-0 decays to other modes, such as to three or five 
pions, which mainly proceed through three gluons. The 
2(7r+7r^) and 3(7r"'"7r~) mass spectra for events of the ISR 
processes e+e" — )■ 2(7r+7r~)7 and e+e" — )■ 3(7r+7r~)7, in 
the mass regions of the J/ip and ip{2S) resonances, are 
shown in Figs. [78] and [79l respectively. From the fits to 
the mass spectra the numbers of J/ip and ip{2S) events 
and also the level of the nonresonant background are de- 
termined. The latter is proportional to the value of the 
nonresonant e+e~ — ?> 2(7r"'"7r~) or e^e~ — » 3(7r^7r~) c ross 
section. In the BABAR paper (jAubert et all l2006b[) the 
ratio 



R 



67r2r( J/V' ^ e+e-)B{J/ip 



f)/m^ 



j/i, 



J/i' 



f("^J/^) 



(19) 



is calculated, where tTe+e-_j./ is the value of the nonres- 
onant cross section to the final state / at the J/i/j mass. 
The numerator of the ratio represents the integral over 
the J/tp excitation curve. The Rj/,p values for the 47r, Gtt, 
2ii:27r, 2i4:47r, and AK final states are hsted in Table HV] 
together with the Rj/^ value obtained for the fJ-^ IJ-^ final 
state. The Rj/,p values for the 47r and Gtt final states are 
closer to that for compared to the final states with 

kaons and indicate that the single-photon exchange dom- 
inates for the J/i/j decays into these modes. For the J/ip 
decays to the final states with kaons, which can contain 



FIG. 79 The 3(7r vr") mass distribution for ISR-produced 
e^e~ — > 3(7 r'^7r~) events i n the m ass region around the J/'tp 
and ip{2S) (jAubert et all l2006bl 'l: there are clear signals at 
the J/-!/) and ip{2S) mass positions. 



a sizeable isoscalar component, the single-photon transi- 
tion is expected to be less dominant, as indicated by the 
larger central values of the ratios. 

TABLE IV Ratios of the J/ip partial production rates to 
contiimum cross sections Rj/^ (see E q.( I19p). The result 
for fi^fi^ is from Ref. (lAubert et all 2004ah. The result 
for 3(7r^7r~) is from Ref. (jAubert et al. I. l2006bl ) and the re- 
sults f or 2(n+n-) K+K- -K+TV- and K+ K' K+ K' are from 
Ref. ( Aubert et aU 12005^ 1. 



Final state 



R. 



2(7r+7r-) 
3(7r+7r-) 
2(7r+7r-7r'') 
K+K'2{n+n~) 

K+K-K+K- 



(MeV) 



85.1 ±7.9 
106 ± 10 

99.1 ±6.5 
122 ± 10 
166 ± 19 
138 ± 32 
84.12 ±2.69 



VI. ISR STUDIES IN THE CHARMONIUM REGION 

In this chapter we will discuss recent progress in 
the charmonium spectroscopy mainly achieved due to 
the ap plication of the ISR method, see also recen t re- 
views ( Brambilla et am201ll : iPakhlova et al. I. I2010D . We 
will start with the description of the open charm final 
states addressing later so-called charmonium-like states, 
presumably states with hidden charm. 
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A. Final states with open charm 

For a quarter of a century our knowledge of the 
vector charmonia above the threshold of open charm 
production (throughout this section referred to as V 
sta tes) was based on th e pioneer expe riments of Mark- 
I (ISiegrist et all Il976[) and DASP (|Brandehk et all 
I1978D . Even such basic parameters of the ip mesons 
as mass, width and leptonic width were known with 
large uncertainties mainly determined by low s tatistics 
of the old experiments. In Ref. (Sethi. |20Q5[ ) an at- 
tempt was made to use the updated informa t ion o n the 
R val ues from Crysta l Ball ( Osterheld et all Il986l ) and 
EES ()Bai et all l2002t ) to improve these parameters. Fi- 
nally, the BES Collaboration performed a global fit of 
the data on R collected by BES in t he energy range from 
3.7 to 5 GeV (|Abhkim et all l2008l ). In some cases the 
obtained values of mass, width and leptonic width for the 

states differ significantly from the older values and still 
suffer from big uncertainties caused by insufficient statis- 
tics and model dependence primarily due to numerous 
thresholds of charm production opening in this energy 
region. It became clear that serious progress would be 
possible after tedious exclusive studies, which recently 
became possible due to ISR analyses of B AB AR and Belle 
based on very large integrated luminosities. 

Exclusive e+e~ cross sections for hadronic final states 
containing charm mesons in the \/s =3.7-5 GeV/c^ energ; 
range were meas ured by BABAR jAubcrt et ai, 2007| 
l2009b[ ) and Belle (jPakhlova et all . [2007, ,2008aiib <ai20& 
using ISR to reach the charmonium region. Note that 
in these analyses Belle systematically employs a partial 
reconstruction technique to increase the detection effi- 
ciency and suppress background. 

The DD cross sections in the entir e charm 
energy range from Belle (iPakhlova et all l2008al ) 
and BABAR (jAubert et all l2007el ) are shown in 
Figs. [80l a).(b) and are consistent with each other. 
Both exhibit clear evidence of structures near 4.1 and 
4.4 GeV/c^. They also observe a structure fFigs. [80la) 
and (b)) at 3900 MeV which must be taken into ac- 
count to describe the DD cross section and R in the 
region between the V(3770) and V(4040). This en- 
hancement is not considered as a new cc resonance, 
as it is qualitatively consistent with the energy de- 
pendence of the sum of the cross sections for various 
channels opening in th is energy range pre dicted in a 



coupled-channel model (lEichten et a/.l. Il980l) . The DD* 
cross sections from B elle ("Pa khlova et al. I I2007D and 
BABAR ( A ubert et aL 2009 b.) shown in Figs. |5gi;c).(d) 
exhibit a single broad peak near threshold (close to 
the ijj j AOAO) position), whe reas the D*D * results from 
Belle (IPakhlova et all iwofl) and BABAR (|Aubert aU 



l2009bl ) ('Figs.[5CTe).(f)) feature several local maxima and 
minima in this energy range. 

BABAR (|Aubert all l2009bD performed unbinned 
maximum likelihood fits to the DD, DD*, and D*D* 
spectra. The expected ip signals were parameter- 



ized by p-wave relativistic Breit-Wigner (RBW) func- 
tions with their parame ters fixed to the PDG08 val- 
ues (|Amsler et all l2008l ). An interference between the 
resonances and the non-resonant contributions was re- 
quired in the fit. The computed ratios of the branching 
fractions for the "0 resonances and the quark model pre- 
dictions are presented in Table Ej The BABAR results 
deviate from some of the theoretical expectations, which 
often differ from each other. 
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FIG. 80 Measured e+e 

sections for y^=3. 7-5.0 Ge V/c^ from B ell e an d 

BABAR, showing (a) DD (IPakhlova et all l2008al') . 

(b) D D (|Aubert ej oZ.l . l2007el ) , (c) D+D" (jPakhlova eTaU . 
l2007l ): (d) DD* for D=D ° (solid squares) and 



.s) 

M 



2009b 



D=D'^ ( open circles) (|Aubert et al , _^ 

(e) D' +P— (IPakhlova eraUl2007l ): (f) D*D * (lAubert et aZ. 
2009bD : (g) D°D-n+ (IPakhlova et al. 



2008bl )^ (h) D°D*- n+ JPakhlova et all |2009| ): 



(i) Ac Ac (|Pakhlova et all 1200801 ). Vertical dashed lines 
indicate 4> masses in the region. 



The 



— » £)°D TT^ cr oss section measured by 

Belle (IPakhlova et all l2008bD is shown in Fig. [SOKg) 
and exhibits an unambiguous ^/'(4415) signal. A study 
of the resonant structure shows clear signals for the 
1)2(2460)° and D2(2460)+ mesons and constructive in- 
terference between the neutral D^D2{2A6Q)^ and the 
charged D^D2{2A60)~^ decay amplitudes. Belle per- 
formed a likelihood fit to the D^D^t:^ mass distribution 
with a ?/'(4415) signal parameterized by an S'-wave RBW 
function. The significance for the signal is ~10(t and the 
peak mass a nd total width a re in good agreement with 
the P pCOe (lYao et all l2006l) values and the BES fit re- 
sults (|Abhkim et all l2008[ l The product of the branch- 
ing fractions S(V'(4415) ^ DD^{2460)) x 6(5^(2460) ^ 
Z?7r+) was found to be between 10% and 20% depend- 



42 



TABLE V Ratios of branching fractions for the ^/'(4040) , 
'0(416 0) and i/)(4415) resonances from BABAR (jAubert et all 
Theoretic a l expe ctation s are from models d enoted 
(|Eichten et all I2006I '). and 



Po ( Barnes et all 



pKp (jSwanson et al. 



20051) 
. I2OO6I ) 



State 


Ratio 


Measured 






pKp 


7/'(404Q) 


DD/DD* 


0.24±0.05±0.12 


0.003 




0.14 




D'D*/DD* 


0.18±0.14±0.03 


1.0 




0.29 


V'(4160) 


DD/D'D* 


0.02±0.03±0.02 


0.46 


0.08 






DD* /D*D' 


0.34±0.14±0.05 


0.011 


0.16 




i/.(4415) 


DD/D*D* 


0.14±0.12±0.03 


0.025 








DD* /D*D' 


0.17±0.25±0.03 


0.14 







ing on the ?/'(4415) parameterization. The non-resonant 
D^D^TT^ branching fraction was found to be <22% of 
e(V'(4415) ^ L>I)*(2460) D°D-n+). Simi- 

larly, the energy dependence of the cross section of 
the D°D*~7r+ final state, shown in Fig. ISOTh) has 
been measured by Belle (jPakhlova et all l2009l ): a 
marginal signal of the ^(4415) is found (S.lcr), and 
its branching fract ion was limited to <1 0. 6%. Very 
recent ly BABAR ddel Amo Sanchez et "o/] . I2OIOI) and 
Belle (jPakhlova et all 120111 ) reported consistent results 
on the cross sections of D+D~, D+D*+ and D*+D*~. 

The Belle collaboration has also me asured the cross 
section of the process e~''e~ — > A+Aj ( Pakhlova et all 
l2008cl) . Because of the large number of the A^ de- 
cay channels with small branching fractions full recon- 
struction of both Ac is not effec tive. The strategy of a 
search for A+A~7 events at Belle (jPakhlova et a/.l . l2008c[ ) 
was the following: one of the Ac baryons was recon- 



structed using three decay modes pKg, pK 



A7r+. 



Then in the spectrum of masses recoiling against the 
A+7 system, a peak at the A^ mass was searched 
for. This peak presumably corresponded to the process 
e+e^ A'^A'^-f. The resulting exclusive cross section 
of the process e+e~ — ^ ^t^c shown in FiglSCTi). 
The cross section is nearly flat from the threshold up 
to 5.4 GeV/c^ except the region just above threshold, 
where a peak with the mass M = 4634;;^° MeV/c2, 
width r = 92^2? MeV/c^ and significance of 8.2 a is 
observed. The state is denoted as X(4630) and the 
product of the branching fractions measured for it is 
B{e+e-) X S(AcAc) = (0.68±0.33) x IQ-^. The nature of 
this enhancement remains unclear. Although both mass 
and width of the A (4630) are consistent within errors 
with those of another Belle state F(4660), that was found 
in il;(2S )TnT decays via ISR and is described in the next 
section / Wang et ai\ . \200l\ ) . this could be coincidence and 
does not exclude other interpretations. 

Although in general the energy behavior of the ex- 
clusive cross sections from BABAR and Belle qualita- 
tively follows the expecta tions of the coupled-channel 
model ( Eichten et al. I. I1980D . some features are not repro- 
duced by theory. This is confirmed by the measurement 



of CLEO (|Cronin-Hennessv et a/.l . l2009f) . which scanned 
the energy range between 3.97 and 4.26 GeV and re- 
ported the cross sections for final states consisting of two 
charm mesons {DD, D*D, D*D*, D+D', D^+D-, a.nd 
D*^D*~) as well as for those in which the charm-meson 
pair is accompanied with a pi on. The updated potentia l 
model predictions of Eichten (|Eichten et all Il980l |2006[ ) 
fail to describe many features of the data. 



B. New charmonium-like states 

The first observation of an unexpected 
vector c harmonium-lik e stat e was made by 
BABAR (|Aubert et all l2005al ) in ISR production 
of y(4 260) — > JMTr+TT" , which was later up- 
dated ([Aubert et all l2008a|) w ith twice the data, 
as sho wn in Fig. [j Tl C LEO (|He et all |2006| ) and 
BeUe (|Yuan et all l2007l) confirmed the BABAR 

broader struc- 
Aside 
the updated 



result, but Belle also found a smaller, 
ture at 4008 MeV/c^, as seen in Fig. 
from the lower mass state, for which 
BABAR ( Aubert et a^.l . [20"08a|) analysis placed an upper 
limit, the three sets of measurements were quite consis- 
tent in mass and width, as shown in Table | VI[ but only 
roughly so in strength. BABAR (jAubert et all l2007al ) 
found one more appare nt enhancement Y (4360) in 
ip{2S)TT+TT- , which Belle ^^ITdl^OO^ measured 
with somewhat larger mass and smaller width, as seen 
in Table I VIII Belle also found a second structure near 
4660 MeV/c^ in the same f i nal st ate, as seen in Fig. [531 
(A combined fit (|Liu et all 120081 ) to Belle and BABAR 
'0(25')7r+7r~ data found consistency between them.) 
Because dipion transitions between vector quarkonia are 
commonplace for charmonium and bottomonium, it was 
natural to ascribe the K's to excited vector charmonia. 
A number of additional features of these states are in 
conflict with this hypothesis. Only one, F(4660), is 
remotely near a predicted 1" cc state (l'^I?i). The 
y(4260) and V(4360) did not show up in inclusive 
hadronic cross section (i?) measurements (jBai et all 
|2002| ). as would be expected of such states (there is no 
fine-grained i?-scan data near y(4660)). 

A comparison of the measured J/^7r"'"7r~ and total 
hadronic cross sections in the ^/s ~ 4260 MeV region 
yields a lower bound for T{Y — >■ J/-!/)7r"'"7r~)>508 keV 
at 90% C.L., an order of magnitude higher than 
expecta tions for conv entional vector charmonium 
states (|Mo et all Hooi). Charmonium would also 
feature dominant open charm decays, exceeding those 
of dipion transitions by a factor expected to be ^100, 
because such is the case for ^(3770) and ip{4160). 
As summarized in Table IVIIIi no such evidence has 
been found, significantly narrowing any window for 
either charmonia or, in so me cases, qua r k-gluo n hybrid 
interpretations. CLEO (jCoan et all l2006( ) studied 
direct production of the y(4260) in e~^e~ collisions and 
identified the only non- J/i/'tt+tt^ decay modes seen 
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FIG. 81 The invariant mass of J/tpTT'^Tr~ candidates pro- 
duced in initial state radiation, e"'"e~ — >■ 7/s_r J/t/jtt+tt" . 
Points with error bars represent data, and the curves 
show the fit (solid) to a_signal plus a linear background 
(dashed) (|Aubert et am2008al ). 




J T-— I— ■i.......,H.».'.5li.bl:iiiiLii.lLtttiL; 



M(7cVJ/^) (GeV/c^) 
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FIG. 82 The invariant mass of J/tpTT'^Tr~ can didates pro- 
duce d in initial state radiation studied by Belle (|Yuan et all 
|2007| ). with J/^/j-sidebands already subtracted, unlike Fig. 1811 
Points with error bars represent data, the solid curve shows 
the best fits to the data to two resonances including inter- 
ference with a floating phase, and the dashed and dashed-dot 
curves show the two pair of individual resonance contributions 
for the two equaUy probable best-fit phases. 



SO far, J/ipTT^Tr'^ and J/tpK^ , occuring at roughly 
half and one-sixth, respectively, of the J/V'tt+tt" rate. 
The J /il}K'^K~ decay mode was also observed by 
Belle (jYuan et a/.l . l2008l) . 

Any interpretation for these vector states will not only 
have to explain their masses, widths, and manifest reluc- 
tance to materialize in open charm or unflavored light 
meson final states. The dipion invariant mass spec- 
tra exhibit curious structures, as seen for the F(4260) 



TABLE VI Measur ed properties of the y(4260) 
J/il)-K^-K~ . The Belle (jWang et ali 120071) single-resonance fit 
result is quoted to allow for comparison to the other two. 



Quantity 



M 

(MeV/c^) 



r 

(MeV) 



B X r„ 

(eV) 



Value 



From (x^/ndf) 



4259±8+6 
4263±6 

-17_ 
-16- 

4263±5 



B ABAR (Aubert et al. 2008a) 
Belle (Yuan et al . 2007) 
CLEO (He et al. 2006) 
Avg (1.8/2) 



88±23t| 
126±18 

73^25 ±5 
108±15 



BABAR (Aubert et al. 2008a) 
Belle (Yua n et al., 2007) 
CLEO (He et al . 2006) 
Avg (2.4/2)^ 



5.5±1.0+o? 

9.7±l.l' 
S.Qtli ± 1.8 

8.0±1.4 



BABAR (Aubert e t al, 2008a) 
Belle (Yuan et al .. 2007) 
CLEO (He et al. 2006) 
Avg (6.1/2) 



TABLE VII Measured properties of the two enhancements 
found in the tp(2S)'K'^ n~ mass d i stribu tion, the 1^(4360) and 
F(4660). Liu et al. (jLiu et ail boOSi ) performed a binned 
maximum likelihood fit to the combined Belle and BABAR 
cross section distributions (Fig. 1831) . 



Quantity 


Value From (x^/d.o.f.) 


M 

(MeV/c^) 


4324±24 BABAR (Aubert et al. 2007a) 
4361±9±9 Belle (Wang et al. 2007) 

4353±15 Avg (1.8/1) 
4355+ 10 ±9 Liu (Liu et al. 2008) 


r 

(MeV) 


172±33 BABAR (Aubert et al. 2007a) 
74±15±10 Belle (Wang et al.. 2007) 

96±42 Avg (6.8/1) 
103+}5±11 Liu (Liu et al. 2008) 


M 

(MeV/c^) 


4664±11±5 Belle (Wang et al. 2007) 
4661+?±6 Liu (Liu et al.. 2008) 


r 

(MeV) 


48±15±3 Belle (Wang et al. 2007) 
42+1^ ±6 Liu (Liu et at.. 2008) 



in Fig. I84l (lAubert et al 



l2008ah . for the F(4360) in 

Fig. ESlJa) (Wang et al!, '20071), and for the y(4660) in 
Fig. imb) (Wang et al, 2007). The first shows a dis- 
tinctly non-phase-space double-h ump structure wh ich is 
qualitatively confirmed by Belle ([Yuan et al ]. l2007l ). the 
second exhibits a majority of events at higher masses, and 
the third indicates a quite dominant /o(980) component. 



VII. SOME IMPLICATIONS FOR THEORY AND 
PERSPECTIVES 

The progress in precision of the low energy data on 
e+e~ — > hadrons achieved recently due to ISR studies 
allows an update of the estimation of the hadronic con- 
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FIG. 83 From a binned maximum likelihood fit (|Liu et all 
120081 ') of combined Belle and BABAR data, the V'(2S)7r+7r" 
invariant-mass cross section as a function of ^/s. The solid 
circles and stars show the Belle and BABAR data, respec- 
tively. The solid curve shows the best fits to the data to two 
resonances including interference with a fioating phase, and 
the dashed curves show the contributions of two pairs of indi- 
vidual resonances for the two equally probable best-fit phases. 
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FIG. 84 The dipion invariant mass distribution on ISR- 
produced y(4260) -i> J/tp-K+n' decays, where points rep- 
resen t data and the line h istogram phase-space MC simula- 
tion (jAubert et aU 1200881 ). 



1 + -K-) 

et all 



TABLE VIII Upper hmits at 90% C.L. on the ratios 
cr(e+e- y T)/a{e+e- -» F -» J/tpTT-^ 

at = 4.26 GeV/c^ (CLEO (|Cronin-Hennessv et 

\200% ) and B{Y T)/BiY J/i >n+n-) (for 

y(426 0)) (BABAR (lAubert et all. l2007el . l2009bD and 
Belle l|Pakhlova et all\200l\ )) . where T is an open charm final 
state. 



DD 
DD* 

D'D* 

DD'n 

D*L»*7r 
DsDs 



y(4260) 



4.0 (Cronin-Hcnncss v et al . 
7.6 (Aubert et al. , 20076) 
45 (Cronin-Hennessv et al . 
34 f Aubert et al. 2009b) 
11 f Cronin-Hennessv et al . 
40 (Aubert et al^ 2009b) 
15 f Cro nin-He nnessv et al. . 
9 (Pakhl ova et al. 2007) 

8.2 (Cronin-Hennessv et al . 

1.3 f Cronin-Hennessv et al . 
0.8 (Cronin-Hennessv et al. 
9.5 (Cronin-Hcnncssv et al. 
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2009) 
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FIG. 85 The dipion invariant mass distribution on ISR- 
produced (a) y(4360) i>{2S)TV+Tv~ and (b) y(4660) 
'4>{2S)n'^n~ , where points represen t data and the lin e his- 
togram phase-space MC simulation (jWang et a^L [20071 ) . 



from a dispersion integral (iBouchiat and Michell . Il96ll : 
iGourdin and de Rafael . Il969l) in which the integrand 
contains a combination of experimental data on cross 
sections of e+e" — >■ hadrons and perturbative QCD. The 
integral ranges from the threshold of hadron production, 
i.e., from the -k^^ threshold to infinity: 



liad,LO 



ds 



R{s) k{s) 



(20) 



tribution to the muon anomalous magnetic moment to 
be performed. It is well known that the precision of the 
Standard Model prediction of this quantity is limited by 
the contributions from strong interactions. These are 
conventionally separated into a theory-dr iven light-by- 
light contribution, see a recent review in ( Prades et all 
|2009[ ). and two experiment-driven vacuum polarization 
contributions, the dominant lowest-order and higher- 
order parts. The lowest-order term can be calculated 



The function K{s) in the integration kernel is rather 
smooth, whereas a factor emphasizes the low-energy 
part of the spectrum. Of particular importance is the 
process e^e~ — )■ 7r''"7r~(7), which provides about 73% of 
the lowest-order hadronic contribution and about 62% of 
its total quadratic error. 

In most cases new ISR results from BABAR are con- 
sistent with previous measurements and have compara- 
ble or better accuracy. However, not always these results 
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Process 


Before BABAR With BABAR 




2.45 ±0.26 


3.25 ± 0.09 


27r+27r" 


14.20 ±0.90 


13.09 ± 0.44 


37r+37r" 


0.10 ±0.10 


0.11 ±0.02 


27r+27r~27r" 


1.42 ±0.30 


0.89 ± 0.09 



TABLE IX The contribution of some multipion processes to 
^had.LO integrated from threshold to 1.8 GeV for the measure- 
ments before BABAR and including the new BABAR results. 
For the it^tv~it^ final state the contribution of the uj and 
mesons is excluded. All values are in units of lO"^''. 



agree with the corresponding old datasets. For exam- 
ple, from Fig. [53] discussed in Ch. 3 it is clear that the 
cross section of the process e+e ~— » 7r+7r~7r° obtained by 
BABAR (A ubert et aU l2004bl) is consistent with that 
of SND (lAchasov et al\ . I2002D below J^_^J^A_GeV, 
but i s much higher than that at DM2 ( Antonelli et all . 
I1992D above this energy. The energy dependence of the 
cross section observed by DM2 is also inconsistent with 
other m easurements (s e e the discussion of this problem 
in Ref. ( Achasov et all [20021) ) and the existence of the 
rather well established a;(1420) and a;(1650) resonances. 
The contribution of this process to a^^'^'^'-', which was 
equal t o (2.45±0.26)-10-i° before BABAR (|Davier et all 
l2003b[ ). becomes (3.25 ± 0.09) ■ 10~^° after the new re- 
sults are taken into account (jDavieil . [2007). For the 
process e^e~— !■ 27r+27r~, which cross section is one 
of the lar gest above 1 GeV, t he new BABAR mea- 
surement ( Aubert et all l2005bl ) is in good agreement 
with the older results and after taking them into ac- 
count the precision of the corresponding contribution 
improves by a factor of two. Another e xample is the 
measu rement of two six-pion final states ( Aubert et all 
l2006bl ). In Figs. |86] and [87] we compare the cross sec- 
tions from BABAR with those in older measurements. 
It is clear that the improvement is dramatical because 
older measurements were too imprecise to make a rea- 
sonable prediction. We summarize the discussed contri- 
butions to a^'^^'^'^ integrated from threshold to 1.8 GeV 
for the measurements before BABAR (see the references 
in Ref. (iDavier all l2003bl) ) and with BABAR in Ta- 
ble [TX] (jDwielliQol . 

The calcul ation using al l multibody niodes nieasured 
by BABAR (l Aubert et al. I, l2004bl l2005bl l2006bl l2007cD 
together with the relevant older measurements (for a 
complete list of references see Ref. (jDavier et aLl . l2003bl )) 
gives for the contribution of hadronic continuum from 
threshold to 1.8 GeV the value (54.2 ± 1.9) • 10-i°. 
It is consistent with the result of the older calcula- 
tion (jDavier et aLl . l2003b[) . (55.0 ± 2.6) • lO"!'', and more 
precise. 

A very important step in the calculations ofa]^'^''"^ has 
been recently made in Ref. ( Davier et aLl . [20l6a[) . which 
for the first time took into account the BABAR mea- 
surement of the reaction e^e" 
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FIG. 86 Th e cross secti on of th e process e" ^e 



27r+27r~27r° (lAubert ail l2006bl : iBacci et all Il981 
Baldini et a lj. Il988l : ICosme et all Il979l : lEsposito et al 



19811 ). 




^ TT+TT (lAubert et all 

l2009al ) and finally in Ref. (jDavier et aZI . l201ll) . which 



FIG. 87 The cross se cti on of t he process e ^e — > ■ 
37r+37r~ (lAubert et all l2006bl : iBaldini et ail 1 19881 : 
iJean-Marie et aZ.l . ll976l ). 



also accounte d for the KLOE re sults on the tt+tt^ fi- 
nal state ([Ambrosino et all 12009.. ^011). It uses the 
whole set of experimental data in the energy range up 
to 1.8 GeV, the region dominated by hadronic reso- 
nances, and perturbative QCD for the contribution of 
the quark continuum beyond that energy. In particular, 
they modify the treatment of the a;(782) and (/)(1020) res- 
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Experiment 


al'"'-^°[nn] [10""'] 


Bcvc [%] 


BABAR 
KLOE 
CMD2 
SND 


514.1 ±3.8 (1.00) 
503.1 ±7.1 (0.97) 
506.6 ±3.9 (0.89) 
505.1 ±6.7 (0.94) 


25.15 ±0.18 ±0.22 (1.00) 
24.56 ± 0.26 ± 0.22 (0.92) 
24.96 ±0.21 ±0.22 (0.96) 
24.82 ± 0.30 ± 0.22 (0.91) 



TABLE X Estimat ed a!;'"''^°[7r7r] and Bcvc contributions 
from the e^e~ data ijPavier et aZ.I . [20111 ). 



onances usins; the data from Ref. ( Aubert et al 
and includ e prelimin a ry da ta of BABAR on e~^e~ — >■ 
7r"'"7r~27r (jPruzhininl . l2007t ). They also use isospm m- 
variance to estimate the contributions from several un- 
measured final states with six pious and KK{m:) relat- 
ing them to those of known chann els. For other channe ls 
they refer to earlier calculations ( Davier et all l2003al lbl: 
lDayieil . [20Q7l) . 

Contributions to a'^'*'*''^'-'f7r7rl from the individual tt+tt" 
cross sections measured at BABAR (lAubert et ali . 

|201lD, 



l2009al ). KLOE (lAmbrosino et all 



(lAkhmetshin et al.. 



CMD2 

SND dAchasov et al. 



2004 



a 



are 



( Davier et 



iven 
201lh 



in the 
. In the 



and 
middle 
right 



column of Table | 
column we list the corresponding CVC predictions for 
the — > TT^vr^i^T- br anching frac t ion co rrected for 



isospin-breaking effects ( Davier et aLl . [2010b[ ). Here the 
first error is experimental and the second estimates the 
uncertainty in the isospin-breaking corrections. For each 
experiment, all available data in the energy range from 
threshold to 1.8 GeV {rrir for Bcvc) are used, and the 
missing part is completed by the combined e+e" data. 

The average of the four sep arate results give s 
^had,LO[^^] = (507.8±2.8tot)-10-i° (iDavier et a/.l.l201ll) 
The co mparison with their previous result ( Davier et al\ . 
l2010bD . al'^^^^^lTTTT] = (503.5 ±3.5tot) ■ 10-^°, shows that 
the inclusion of the new BABAR data significantly in- 
creases the central value of the integral. For the higher- 
order hadronic contributions to a'^'*^''"'-*, which are also 
estimated based on e^e^data, there is a slight gain in 
accuracy from — 9.79±0.08cxp±0.03rad (^Hagiwara et aZj, 
l2007l ) to -9.79±0.06cxp±0.03,.ad (Hagiwara et al, 201l|. 

A compilation of recent results for a^, from 
which the central value of the experimental aver- 
age ( Bennett et [2006) has been subtracted, is 
given in Fig. [551 The shaded vertical band indi- 
cates the experimental err or. The SM prediction s are 
tak en from: HMNT 07 (iHaeiwara et all |2007[ ). JN 
09 (I Jeeerlehner and Nvffcler'. '20091) . Davier et al. 09/1 
(r-based) (|Davier et al . .2010b.) , Davier et al. 09/2 



(e+e~ -based, with BABAR 
l2010al ). HLMNT 11 



data) (|Davier et all . 
-based, with BABAR and 



KLOE TT+TT" d ata) ([Hagiwara et all [201 1[ ). DHMZ 10 
([Davier et aL[ [201ll). BDDJ 11 ((r and 



(r and 



) (|Benavoun et all\201lh 



There have been only a few tests of CVC based on the 
new data. The most interesting final state is of course 
that with two pious, where a serious discrepancy be- 



HMNT 07 (e*e ) 
-285 + 51 

JN 09 (eV) 
-299 ± 65 

Davier etal. 09/1 (x) 
-157 ±52 

Davier et al. 09/2 (e*e" ) 
-255 ± 49 

DHMZ 10 (T) 
-195+54 

DHMZ 10 (eV) 
-287 ± 49 

HLMNT 1 1 (e*e") 
-261 ±49 

BDDJ 1 1 (e*e + t) 
-335 ±53 

BNL-E821 
0±63 



-800 



-600 



-400 



-200 
(a -a„ )xlO 



11 



FIG. 88 Compilation of recent results for a^. 



tween the e+e and r data was reported ( Davier et all 
l2003al[bl ). The results of the most recent tests in this 



channel ( Davier et al I [20I0bl [20 III) are based on the 
reevaluation of isospin-breaking corrections. The pre- 
dictions for the branching fraction of — >■ Ti^iT^Vr 
shown in Table [X| can be co mpared to the world av - 
erage value of (25.51 ± 0.09)% (|Nakamura et all [2OI0l ). 
Although the difference between the CVC prediction 
and the experimental value i s less significan t than pre- 
viously (|Davier et all [2003b[ : [Davier[ . [20071) . it is still 
substantial for all groups but BABAR. A new ap- 
proac h to the problem was s u ggest ed very recently in 
Ref. ( Jegerlehner and Szafrod . 20Ilh . where the p — 7 
mixing was properly taken into account. The CVC 
prediction for the branching fraction to two pious is 
(25.20 ± 0.17 ± 0.28)% in good agreement with t he di- 
rect measurement. Finally, Ref. ( Benavoun et aL[ . [201lh 



reports consistent results on e^e~ and r in the 27r chan- 
nel and obtains a theoretical prediction for a^, which is 
lower than the experimental value by 4.1cr. 

One more test of CVC that included recent 
ISR results from BABAR ha s been performed in 
Ref. ( Cherepanov and Eidelman[ . [20091) . The authors 
use CVC together with the data on e+e^ — > t^tt+tt^ 
and e'^e^ — > r]'TT~^TT~ to estimate the branching 
fraction of the corresponding r decays. For the 
former final state the e stima t e based on the older 
data d Akhmetshin et al\, '2000; 'AntoneUi et dl, [I988I : 
[Delcourt et al.. X982; Druzhinin et al., 1986) predicts 
for the branching fraction (0.132 ± 0.016)%, somewhat 
smaller but not incompatible wit h (0.165 ± 0.0 1 5)% ob - 
tained from the BABAR data (|Aubert et all [2007d ). 
The average of the two gives the CVC prediction of 
(0.150 ± 0.016)%, in good agreement with the new world 



47 



average B(t- rjir' Vr) = (0.139 ± 0.0101% that 
uses the new presise measurement at Belle ( Inami et all 
l2009f ). For the B{t- Ty'Tr^TrOj/^) they give an up- 
per limit of < 3.2 • 10"^ at 90% CL, which is a fac- 
tor of 2.5 more restrictive than the upper limit based 
on th e only existing mea surement by CLEO: < 8 
10-5 (iBergfeld et aLl . ll997D . 



There are two recent evaluations of Aa|j^j(Mf ), the 
hadro nic contribut i on to the running a f r om fi ve fla- 



vors 



fPavier et al 



Ref. (Hagiwara et al 



201lt iHagiwara et all l201l[ ). In 

201 1[ ) a data set of e"'"e~cross sec- 



tions includes multibody data from BABAR and 27r 
data from KLOE and the calculation gives the value 
0.02760 ± 0.00015, slightly higher and significantly more 
accurate than the previously accepted value 0.02758 ± 
0.00035 (iBurkhardt and Pietrzvfll2005l) . The estimation 



performed in Ref. (jPavier et al. 



20 111) additionally uses 



the BABAR data on the tttt final state and perturbative 
QCD between 1.8 and 3.7 GeV and gives an even more 
precise value 0.02749 ± 0.00010. 



VIII. CONCLUSIONS 



Only part of the available ISR data sample has been 
processed, e.g., for BABAR it is about 1/3, analysis is 
in progress. Belle has only started a corresponding data 
processing. 

If existing projects of Super B-Factories are approved, 
prospects of reaching an integrated luminosity by a factor 
of 30-100 exceeding that today appear. Such experiments 
will improve accuracy for many processes which studies 
are now statistically limited. 
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Successful experiments at high-luminosity e^e^ collid- 
ers (0- and -B-factories) opened a new era in a study of 
e"'"e~ annihilation into hadrons at low energies using a 
novel method of initial-state radiation usually referred 
to as ISR or radiative return. 

Modern detectors operating at these factories which 
collected unprecedentally high integrated luminosity al- 
low this method to compete with direct e'^e^ experi- 
ments. 

A lot of new data on the cross sections of e+e^ an- 
nihilation into hadrons were obtained using ISR, first 
of all, on the production of mesons from threshold of 
their production ^ 2r7i7r to the cm. energy of about 4-5 
GeV. More than 30 processes have been studied in which 
mesons and hadronic resonances were produced, many of 
them for the first time. 

Valuable information on the particles with mass of 
about a few GeV has been obtained, primarily on excited 
vector mesons, radial and/or orbital excitations. Param- 
eters of vector charmonia were investigated, new data on 
more than 40 decay channels were obtained, many decays 
observed for the first time. 

New data on production cross sections were obtained 
for various baryons: proton. A, and Ac hypcron, open- 
ing new possibilities for testing form factor models. 

New states (p(1900), y(2175), r(4260), r(4320) . . .), 
some of them with presumably exotic quark structure, 
have been discovered. Their nature is not yet established 
and widely discussed. 

New values of the cross sections obtained using ISR can 
be used for more precise predictions of the muon anoma- 
lous magnetic moment, running fine- structure constant 
at the Z boson mass, tests of CVC and many other the- 
oretical models. 
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